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O ■ Abstract 
(N , 

I show how prior work with R. Wald on geodesic motion in general relativity can be generalized 
to classical field theories of a metric and other tensor fields on four-dimensional spacetime that 1) 
ly-^ ■ are second-order and 2) follow from a diffeomorphism-covariant Lagrangian. The approach is to 

consider a one-parameter-family of solutions to the field equations satisfying certain assumptions 
designed to refiect the existence of a body whose size, mass, and various charges are simultaneously 



, scaled to zero. (That such solutions exist places a further restriction on the class of theories to 

. which our results apply.) Assumptions are made only on the spacetime region outside of the body, 

> 

lO ' so that the results apply independent of the body's composition (and, e.g., black holes are allowed). 

' The worldline "left behind" by the shrinking, disappearing body is interpreted as its lowest-order 

in 

^Sj . motion. An equation for this worldline follows from the "Bianchi identity" for the theory, without 

o . 

, use of any properties of the field equations beyond their being second-order. The form of the 
force law for a theory therefore depends only on the ranks of its various tensor fields; the detailed 



^ ■ properties of the field equations are relevant only for determining the charges for a particular body 

H : 

(which are the "monopoles" of its exterior fields in a suitable limiting sense). I explicitly derive the 
force law (and mass-evolution law) in the case of scalar and vector fields, and give the recipe in the 
higher-rank case. Note that the vector force law is quite complicated, simplifying to the Lorentz 
force law only in the presence of the Maxwell gauge symmetry. Example applications of the results 
are the motion of "chameleon" bodies beyond the Newtonian limit, and the motion of bodies in 
(classical) non-Abelian gauge theory. I also make some comments on the role that scaling plays in 
the appearance of universality in the motion of bodies. 
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I. INTRODUCTION 



In special relativity, a non-interacting body moves in a straight line. Therefore, it is not 
surprising that in general relativity an "infinitesimal test body" (i.e., a body small enough 
that the curvature of the external universe can be neglected, and weakly-gravitating enough 
that curvature it generates can be neglected) will move locally in a straight line, i.e., it 
will follow a geodesic. But from this perspective it does seem quite surprising that strong- 
field bodies like neutron stars and black holes in fact also move on geodesies (in the limit 
of small size). After all, no matter how small or light such a body, the local spacetime 
metric will differ significantly from that of fiat spacetime, and one would therefore expect 
that nonlinear gravitational dynamics — certainly not special relativity — would principally 
determine its motion. Furthermore, since the metrics of different strong-field bodies will 
differ greatly from each other, one would perhaps expect there to be no universal law for 
the motion of strong-field bodies at all. Indeed, the natural assumption would seem to be 
that the motion of a strong-field body depends in detail upon its composition. 

This expectation is incorrect for a very counter-intuitive reason: in general relativity, 
the motion of a small body is in fact completely determined by field dynamics outside of 
the body. This surprising fact was first demonstrated by Einstein, Infeld and Hoffman 
and has become the foundation of a more modern approach to motion termed "matched 
asymptotic expansions" ^ (see also js-S]). The basic physical requirement of this line of 
work is the existence of a region (the "buffer zone") sufficiently far from the body that the 
body field may be approximated as a multipole series, yet sufficiently close to the body 
that the field of the external universe may be approximated in an ordinary Taylor series. 
The vacuum gravitational dynamics taking place in this region then suffice to determine the 
motion. 

A primary purpose of this paper is to determine to what extent this conclusion general- 
izes to other classical field theories. To investigate this question I generalize the approach 
taken in j4| to deriving geodesic motion in general relativity^] In the formalism of 4| a 
small body is characterized by a one-parameter-family of solutions to the vacuum Einstein 
equation describing the region outside of a body that shrinks to zero size and mass with the 
perturbation parameter, A. A family with such behavior is considered by demanding the 



-'^ I do not treat self- force corrections, which were the primary focus of [J]. 
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existence of a second, "scaled" limit wherein the coordinates and metric are rescaled such 
the body is held at fixed size and mass. At A = in the original limit the body disappears, 
leaving behind a smooth spacetime with a preferred worldline, 7, picked out; this worldline 
is interpreted as the lowest-order perturbative motion of the body. We showed that 7 must 
be a geodesic by applying the Bianchi identity to an effective point particle description that 
(remarkably) emerges at first order in A. In this paper I generalize the approach to theories 
that 1) follow from a diffeomorphism-covariant Lagrangian, ensuring a "Bianchi identity" 
and 2) have second-order field equations. 



For this class of theories the method of {4] gives an equation for 7 that depends only 
on "buffer zone" field properties, showing that the Einstein-Infeld-Hoffman idea remains 
correct in a more general context. More specifically, the equation involves, in addition to 
the value and first- derivative of the external fields at the location of the body, various charges 
(understood to include mass as the charge associated with the metric) that are determined 
from the body's fields in the scaled limit (they are "field monopoles"). The results rely only 
on properties 1) and 2) above and are therefore surprisingly independent of the details of the 
theory. In particular, the force law depends only on the form of the Bianchi identity, which 
in turn depends only on the ranks of the tensor fields considered (although extra identities 
following from gauge symmetry can greatly simplify the results). Therefore, the expression 
for the force in terms of the charges and external field values is in fact identical across 
theories with the same types of tensor fields and the same gauge symmetries. However, the 
charges associated with a particular body composition will differ in different theories, since 
the relationship between a given source and the field monopoles it generates will depend on 
the field equations. In this sense — and only in this sense — are there "differences in motion" 
among theories of our class that have the same types of tensor fields and the same gauge 
symmetry. 

In interpreting the results it is useful to distinguish varying degrees of "universality" in the 
motion of small bodies. In the case of general relativity, all small bodies move on geodesies, 
so that their internal structure is completely irrelevant to their motion. In Einstein-Maxwell 
theory, a single number characterizing the body (the charge-to-mass ratio) determines how it 
will move, so that the internal structure is minimally relevant. In scalar-tensor theory, a free 
function of time (the charge-to-mass ratio of the non-conserved charge) specifies the motion 
of a body, so that the internal structure is somewhat relevant. In higher-rank theories a 



finite number of free functions of time characterize the motion of a body. Of these results 
only geodesic motion in general relativity is truly universal in that it applies to all bodies; 
however, I will refer to all of the above results as "universal behavior in motion" , since the 
information required to determine the motion of a small body is reduced from the complete 
description of the body to the knowledge of a finite number of parameters at each time. To 
adopt the language of condensed matter physics, there are thus large "universality classes" 
of small bodies that move in the same way. 

The content of this paper is as follows. In section [TTl I summarize the formalism of {4] 
to derive geodesic motion in general relativity. In section IIIII I generalize the formalism to 
Einstein-scalar and then more general scalar-tensor theories, deriving the scalar force law. 
Note that mass evolution always occurs, and the scalar charge evolution is unconstrained. I 
discuss the results in the context of specific scalar-tensor theories and comment on scaling 
and universality. In section IIVI I apply the formalism to vector-tensor theories to derive the 
vector force law. This surprisingly complicated equation simplifies to the Lorentz force law 
in theories with the Maxwell gauge symmetry. I also derive the simplified force law in the 
case of non-Abelian gauge theory. Finally in section |V] I give the proof that universality in 
motion is achieved via buffer zone dynamics for tensor fields of arbitrary rank. A definition 
and disambiguation of scale-invariance is given in an appendix. 

I use the conventions of Wald and work in units where G = c = 1. Early- alphabet Latin 
indices a,b,... are abstract spacetime indices, while Greek indices /i, i/, ... give tensor com- 
ponents in a coordinate system. When working in coordinates {t,x^,x'^,x^), mid-alphabet 
Latin indices i,j, ... denote spatial components 1 — 3, while a zero denotes the time com- 
ponent t. Mid-alphabet capital Latin indices I, J, ... label members of a collection of tensor 
fields. 



II. REVIEW OF FORMALISM: MOTION IN GENERAL RELATIVITY 



In this section I review the derivation of geodesic motion given in ^ . While the treatment 
here is self-contained, the reader is referred to j4| for more details and significantly more 
motivation. Note that many of the arguments given here will hold identically or analogously 
for the more general theories treated in later sections, in which case those arguments will 
not be repeated. 



The basic approach to motion is to formahze the notion of a "small body" by considering 
a one-parameter-family of solutions to Einstein's equation that contains a body that shrinks 
to zero size with the parameter A. While no universal behavior in motion (nor even any 
definition of "position"G) can possibly be obtained at any finite A, in the limit A — ?■ one can 
hope for a simplified description, whose observables will then approximate observables at 
small but finite A fl The task is therefore to develop assumptions on a one-parameter-family 
to the effect that it contains a body shrinking to zero size. The first realization is that the 
body must also shrink to zero mass, since (roughly speaking) no body can be smaller than 
its Schwarzschild radius. The body will thus disappear in the limit, but it will leave behind 
a preferred worldline, 7, characterizing its motion. Our method of considering such a body 
is essentially to demand that if we zoom in on the presumed shrinking body, then a body is 
recovered. This zooming process is accomplished via the notion of a scaled limit, defined as 
follows. Consider a one-parameter-family of metrics gabW, whose metric components g^ui'^) 
are given in some particular coordinates (t, x*). Introduce both a rescaled metric (jabW = 
^^^QabW and, for a particular time to, rescaled coordinates (t,a;*) = {{t — to) /\,x^ / X). 
Then, the scaled limit is given by 

^JS(^o) = lini^^p(A;to), (1) 

A — ¥\j 

where the limit is taken at fixed scaled coordinate. In this notation, the bar on the "(yf" 
indicates that the rescaled metric is being considered, while the bars on the coordinate com- 
ponent indices "/x" and "z/" indicate that the components of Qab in the rescaled coordinates 
are being considered. (I will continue to adopt this notation throughout the paperO) 

This limit has the interpretating of "zooming in" because a fixed-x'^ observer moves ever 
closer to the shrinking body while the rescaled metric keeps distances finite. A simple 
example to keep in mind is the family of Schwarzschild deSitter metrics of mass parameter A 
^, which clearly contains a shrinking body of the sort we want to consider. In the ordinary 
limit the body disappears, leaving behind the background spacetime of deSitter. But in 
the scaled limit the background "disappears", leaving behind the Schwarzschild metric for 

^ Consider, for example, the impossible task of assigning a center of mass position to a black hole. 
Of course, the value of A at which the physical spacetime is embedded into the one-parameter-family is 

arbitrary. What matters for the application of the simplified description is that corrections to the relevant 

observables are numerically small. 
^ Note also that coordinate indices will always refer to the original cartesian-like (t, x'), even if the coordinate 

componants are being viewed as functions of other variables, such as spherical coordinates. 
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all to- For all one-parameter-families in this paper, I will refer to g'^J = g^u{\ = 0) as 
the "background metric" and to gp^p{to) as the "body exterior metric" (no interior is ever 
considered). The existence of both original and scaled limits reflects the presence of a body 
with the appropriate scaling, and will be our assumptions i) and ii), below. 

An additional assumption is required. To arrive at this assumption, note that the rescaling 
of the metric by effectively cancels powers of that arise in changing to the rescaled 
coordinates, so that one has the simple formula 

^/ip(A; to; t, x') = g^^{\] t = tQ + At, x' = Ax*), (2) 

relating barred components of the barred metric, ^^p, to corresponding unbarred coordinates 
of the unbarred metric, (7^1^. That is, one simply "plugs in" t = to + At and = Ax* to 
compute (jfijy from g/j^j,. This formula shows that if we regard g^^ as a function of new variables 
a = r,/3 = r/X,d,(j) (at fixed t, and with r, 6, (p defined relative to x* in the usual way), then 
the scaled limit ([1]) is given by the limit a — )■ at fixed /3 of the original metric components 
g^^ij. Similarly, the original limit is given by the limit /3 — > at fixed a. Demanding that 
both limits exist is thus the statement of separate continuity in a and /3. A natural extension 
(argued for at length in \^ on the grounds that it excludes certain pathological behavior) 
is to demand joint continutity and in fact joint smoothness in a and /3 (although only 
is required here). This will be our assumption iii), below. Note that the electromagnetic 
analog of iii) has been shown to hold for a family of shrinking charge-current sources in flat 
spacetime 9|]. 

The assumptions of this section for the metric family gab{X) are that there exists coordi- 
nates (t, X*) at each A > such that the following hold. 

• i) For r > A^ for some constant R, the metric components g^v{\) satisfy the vacuum 
Einstein equation and are smooth functions of (t,x*. A). The worldline, 7, defined by 
A = X* = is timelike. 

• ii) The scaled metric components g/iuito] ^) are smooth functions of (A, t, x*) for f > R. 

• iii) The metric components g^i, are smooth functions of (a, /3) at (0, 0) for fixed (t, 9, 0). 

Assumption i) establishes our domain r > A-R and provides the requisite smoothness for 
perturbation theory on that domain. It also lays the groundwork for the interpretation of 



the domain as the exterior of a shrinking body by taking 7 to be timehke. Assumption 
ii) estabhshes this interpretation according to the ideas of the scaled hmit, and assumption 



iii) adds additional "uniformity" properties {4]. Note, however, that this latter assumption 
has an important physical consequence (i.e., it places an important restriction on the type 
of spacetime for which our approximate results will be useful). It requires that there be a 
spatial region both far enough from the body that its field can be approximated in a series 
in inverse powers of distance (corrections in /3 near zero), and close enough to the body 
that the field of the external universe can be approximated as a series in positive powers of 
distance (corrections in a near zero). Therefore by seeking one-parameter-families containing 
a shrinking body, we in fact end up with a mathematically precise version of the usual "buffer 
zone" assumption of the Einstein-Infeld-Hoffman approach and its descendants. No further 
assumptions beyond i), ii), iii) (and their analogs for other theories) are made in this paper. 

Smoothness in a and /3 allows us to taylor expand in these variables to any finite order. 
However, to derive geodesic motion we in fact require only a single derivative in /3, 

g^^{\ t, r, 0) = fe^,(t, 0) + c^,(t, 0)/3 + 0{a) + ©(/J^) 

= 6^,(t, B, 0) + c^,(t, ^, 0)^ + 0(r) + O (3) 

where Oir) near zero is at fixed A/r, and 0{\/r) near zero is at fixed r. Sorting into powers 
of A and r, we have 

g^,{\, t, r, B, 0) = b,,{t, B, 0) + 0(r) + A (c^,{t, ^, 0)^ + Oir')\ + 0{\^), r > (4) 



where the order symbols are for small r and A. From this expression it is easy to read off 

7/11/ 1 A=0; 



series expressions for the background metric gj^J and linear perturbations g'^^J = d\g^ 



gf}{t,x') = h,,{t) + 0{r\ r>0 (5) 

= c,At, OA)l + 0(1), r > 0. (6) 

Since gjj^J is assumed smooth everywhere, b^i, cannot depend on angles and we have written 
h^yif). (This quantity is usually taken to be rj^y by coordinate choice.) We are also interested 
in the consequences of equation ([3]) in the scaled limit. Using equation ([2]), we have 

^^p(A;to;t,r,^,0) = 6^^(to + At,^,0) + c^^(to + At,^,0)^ + O(Af) + O f^j , (7) 



so that the hmit A — )■ at fixed (t, x*) gives 

gtkio] x') = b^,{to) + c^,{to, O,(p)l + (J^^ , (8) 

and the body exterior metric (jp^p is seen to be stationary (independent of i) and asymptot- 
ically flat (constant as f — )■ oo). (Stationarity follows from smoothness of g^ui^) in ^0 This 
supports the idea that gj^p characterizes the body exterior as it would appear in isolation at 
time to- Equations ([5]), ([6]), and ([8]) are key consequences of our assumptions. 

Geodesic motion is now derived as follows. By assumption i), g^uW is smooth in A at 
A = for r > 0. Thus ^f^!] satisfies the vacuum linearized Einstein equation about g''^\ 

G2[<7«] = 0, r>0. (9) 

Now regard distribution on defined on the background spacetime including at 

r = 0, which is possible because its "most singular" behavior is only 1/r. Since IS a 
second-order, linear partial differential operator, it follows from equations ([6]) and ([9]) via 
the analysis of appendix |A] that, distributionally, we have 

Gl%^'^] = N,^it)S^'\x^), (10) 

for some Nf^,y defined on the worldline = 0. (This result is analogous to the well-known 
fact that V^(l/r) = —4:Tt6^{x^). If the explicit form of G^f^il is used, a formula for may be 
obtained for A^"^^ in terms of angle averages of c^^ and its first angular derivatives.) Thus, 
an effective distributional stress-energy of 1/8ttN^^6'^{x^) has emerged at first-order in per- 
turbation theory, supported on the worldline 7. This is remarkable, given that any true 
stress-energy associated with the body is confined to r < A-R and excluded from considera- 
tion; and furthermore, the body need not be "made" of stress-energy at all (as in the case 
of a black hole). 

The strategy is now to apply "conservation" to the "stress-energy" . That is, because the 
linearized Bianchi identity V^G^),^ [(7^^-*] = holds as an identity on all sufficiently smooth 
9ab (^°^ necessarily satisfying the linearized Einstein equation), the distributional linearized 
Bianchi identity also holds as an identity on distributional g^^J , and we must have 

{%At)5^'\x')) = (11) 

in the distributional sense. Here Vq is the derivative operator associated with the background 
metric g'^fj. The consequences of this equation can be determined in a variety of ways. I 
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will proceed by adopting the specific coordinate choice of Fermi normal coordinates (see, 
e.g., jsl) for the background metric gj^J . On the worldline = 0, the metric components 
are Minkowski {g^,^ = t]^^), and the Christoffel symbols are given by Fqq = Fq^ = aj, where 
Oj are the spatial components of the four-acceleration of the worldline (the time component 
Oq is zero). We then have 

{N^o5^^\x)) = S^'\x) [-doNoo + a,iV,o] + d,S^'\x) [iV,o] (12) 
(iV^.(5(3) (f )) = 6^'^ (x) [-doNo^ + a,N.,, + a.iVoo] + S.^^^^ (f ) [iV,,] , (13) 

with repeated spatial indices summed. The coefficients of 5^'^\x) and di5^'^\x) must sep- 
arately vanish, giving NiQ = Nij = 0, as well as ^oA^^oo = and ajA^oo = 0, i.e., geodesic 
motion when A^^oo is non-zero. 

We can interpret A^^oo through its appearance in equation ffTOj) . Since it multiplies the delta 
function, it is clear that A^oo will determine the singular behavior of the metric perturbation, 
i.e., it will determine the coefficient c^j, in equation Furthermore, if one expands the 
background metric according to ([S]), it is clear that only the constant term b^i, (which here 
equals r^^j, by coordinate choice) is relevant for the determination of c^^ via (fTOj) . Thus, 
we may compute this coefficient by using the stationary linearized Einstein equation off of 
fiat spacetime in global inertial coordinates. The solutions are well known and one obtains 
for the time-time component (which is all we need) that 47rcoo = A"oo (so that cqo is in fact 
independent of angles in these coordinates). Observing the appearance of c^i, in the series 
for the body exterior metric ([H]), we conclude that A'qo is 8tt times the ADM mass of the 
body exterior metric. Therefore we define M = {1/8tt)Noq and refer to this quantity as 
the mass of the body. This explains the role of the requirement that A'oo 7^ for geodesic 
motion to hold: there must actually be a body present in the one-parameter-family for the 
curve to be necessarily geodesic. 

Equation (fTOj) for the effective stress-energy may be clarified by introducing M and by 
rewriting in covariant form. Since = (1,0) in Fermi normal coordinates, we have from 
Nij = NiQ = that Nab = MuaUf,- The spatial coordinate delta function becomes a worldline 
integral of the "invariant" four-dimensional delta function ^4(0:, x') = ^il^^£L!!)_ Thus we 
have 

G^ailig^''^] = SttM [ UauM^, 4r))dr, (14) 



where the mass M is constant. This equation shows that the metric perturbations for our 
family are in fact sourced by the usual "point particle" stress-energy (see, e.g., 3|). Thus, 
despite the fact that point particles do not make sense in general relativity j?!, we have 
shown that they do emerge as part of a (mathematically rigorous) approximate description 
of the metric of an arbitrary small body. Furthermore, the "particle mass" M is indeed the 
ADM mass of the body (as measured in the scaled limit). 

The results of this section (i.e., the results of sec. IV of j4|) may be summarized as 
follows. Consider a one-parameter-family of spacetimes containing a body whose size and 
mass decrease to zero, according to the stated assumptions. Then, the ADM mass M of the 
body exterior metric is a constant independent of time to, and, if M 7^ 0, the worldline 7 left 
behind after the body disappears is a geodesic of the spacetime g'^^^ left behind. Furthermore, 
the far-field effective description in linearized gravity is that of a point particle of mass M. 
These results show, in essense, that small bodies move on geodesies while keeping their ADM 
mass constant and sourcing linear perturbations reflecting a point particle of that mass. 



III. SCALAR-TENSOR THEORIES 

A simple generalization of general relativity is the addition of a scalar field. I will first 
consider the ordinary Einstein-scalar theory in detail. I will then discuss the general case, 
which in fact follows from the computations already done. Finally I make some comments 
on scaling and universality. In this in later sections it will be convenient to use a Lagrangian 
formulation. I will use the definitions and conventions of appendix E of Wald [6], except 
that I will denote his fixed volume element e by "d^x" . 



A. Einstein-scalar theory 



The action for general relativity plus a minimally-coupled massless scalar field is given 



by 



S = j (fx^ {R - 2g'''Va4>Vb4>) , 



(15) 



where R is the Ricci scalar constructed from gab- I have chosen the relative normalization so 



that the theory reduces to that of Quinn 



in the appropriate limit. It is helpful to define 
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E^^^l = {-g)-^/^5S/5g''^ and E^^^ = {-g)-^/^5S / 5<p, which evaluate to 

= Ga,-2 (^V,</)V,</. - \gaW"'^c<pVd<P^ (16) 

= Ag'^'VaVkcf). (17) 

The equations of motion for the Einstein-scalar theory are then simply E^^l = and E^'^'^ = 0. 

The formalism requires a "Bianchi identity" for this theory. To derive such an identity, 
consider the variation of the action (1151) with respect to an infinitessimal diffeomorphism. 
Since the action is diffeomorphism- invariant, the variation must vanish, and one has 







J d^x^g {^^(-2^^^') + E^^^-Vact] , (18) 



where is an arbitrary vector field. Integrating the first term by parts, we derive 



S = -^i^'^lV,0. (19) 



VE 

2 

The field equations were not used in deriving this equation, which therefore holds as an iden- 
tity on all sufficiently smooth {gab-, 0}- (This can also be easily checked by direct calculation 
using equations ( fT6l) and (|T71) .) To interpret this identity, note that nonzero values of E^^^ 
and E'^ would normally be interpreted as stress-energy and scalar charge density (respec- 
tively) associated with some matter field. This equation gives the precise non-conservation 
the matter stress-energy in terms of the matter scalar charge density necessary for consis- 
tent coupling of that matter to the Einstein-scalar theory. Although we will always impose 
E^^l = and E"^ = at finite A, non-zero values will emerge in the linearized, distribu- 
tional description (analogously to equation (fTOj) in general relativity), reflecting an effective 
stress-energy and scalar charge of the body. 

We now seek to generalize the assumptions used in general relativity to the Einstein- 
scalar theory. The main requirement is to take the scalar charge to zero along with the size 
and the mass, in order to keep the energy in the field finite. Thus we seek a one-paramater 
family with scalar field behavior like ~ A/r. To characterize this by the existence of a 
scaled limit, the appropriate rescaling (after changing to scaled coordinates) is simply (p = (p. 
(That is, no rescaling is required; however, we still define for notational consistency.) Then 
for an arbitrary family we define the scaled limit as in ([1]), 

E lim0(A), (20) 
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where the hmit is taken at fixed scaled coordinate. I will refer to cj)^'^^ as the body exterior 
scalar field in analogy with the body exterior metric g^p. 

One can now follow the same path of reasoning as in section [Tll leading one to assume 
the existence of original (i) and scaled (ii) limits, as well as the uniformity condition (iii). In 
other words, the appropriate assumptions for this section are those of HIl where the metric 
(?afe(A) is replaced by the pair {gab{\), 4>{\)}, and the required equations are not Einstein's 
equation but the Einstein-scalar equations e]^^ = and E^^^ = 0. (The same coordinates 
(t,x*) and hence worldline 7 are used for the metric and scalar field.) Note that the new 
assumption iii) will imply the existence of a "buffer zone" for the scalar field as well as the 
metric. 

The steps of the derivation of motion now follow completely analogously. The analog of 
(|3]) holds for the scalar field, which leads to the analogs of ([5]), ([6]) and (|8]), given by 

=6M(t) + 0(r) (21) 
0(1) =c[<^](t,^,0)l + O(l) (22) 
0(0) = + cM(to, 9,<p)^ + (23) 

for some fe'"^' and c^'^l (Here 0(^) = 9a0|a=o is the scalar field perturbation, so that one has 
(/)(A) = (/)(°) + Xcp^^^ + O(A^) for r > 0.) The body exterior scalar field (^C) is seen to be 
stationary and to have smooth behavior in 1/f as f — ?► 00, supporting the name we have 
given it. Note, however, that in constrast to the case of the metric, the body exterior scalar 
field retains some "memory" of the external universe, since b^'^^ contains physical information 
about 0(°). As before, the perturbations 0(^) and gjn) satisfy the linearized field equations 
at r > 0, 

E^P[g^'\<^^'^] = 0, r>0 (24) 

^[<^](i)[^(i)^0(i)] = 0, r>0 (25) 

where eI^^^^^ and £^['?^l(i) denote the linearizations of eI^^ and E^'^\ respectively, off of the 
background {g^'^\ 0(o)}. Since the field equations are second order in gat and in 0, the linear 
operators E^^^^^^ and i5['^](i) are also second order in g^^^j and (f)^^\ Thus, as in the case of pure 
gravity, the "1/r behavior" of equations ([6]) and (!22l) implies via the analysis of appendix 1X1 
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that, distributionally, we have 

Elf'V\<P^'^] = Nlf:iit)5^'\x^) (26) 
^M(i)[^(i)^0(i)] = ArM(t)5(3)(3,^)^ (27) 

for some Njf} and N^'^\ As discussed above, a non-zero value of i?|fj would normally be 
interpreted as a matter stress-energy source, whereas a non-zero value of E^'^'^ would be 
interpreted as a scalar charge density source. Thus, despite the fact that no matter stress- 
energy or scalar charge density was considered, distributional descriptions of both have arisen 
effectively in perturbation theory. We can take advantage of this remarkable occurence, as 
before, by using the "Bianchi identity" for the theory. That is, since the linearization of 
equation ([19]), V^jfj^^^ = -l/2E^'f'^^^'^V^(f)^^\ holds as an identity on all sufficiently smooth 
{g'^i'', 0^^-'}, it must hold as an identity on distributions as well. Thus we must have 



(28) 



in the distributional sense. Here Va is the derivative operator associated with the background 



metric g^^fj . Adopting Fermi normal coordinates as before, this becomes 








6^^\x) 



1 



+ d,6^^\x) 



iO 



+ di6^''\x) 



[g] 



from which we determine N} 



[g] 



= as well as 



(29) 
(30) 

(31) 
(32) 



We can interpret A^^qq (t) and N^'^^it) as follows. The formula for E^^, equation (fT6!) . demon- 
strates that the linearization, e]^^^^^ will depend on second-derivatives of g^^fj only through 
the linearized Einstein tensor G^^ij , and will contain no second-derivatives of 0*^^-' . Therefore, 



the identical arguments from the case of pure gravity carry over, and we have (l/47r)A", 



Coo = 2M, where M is the masj^ of the body exterior metric at time to = t. Similarly, from 



[g] 

00 



^ Applying the usual notion of mass to Einstein-scalar theory makes sense because the scalar field stress- 
energy is quadratic in first-derivatives of (f), so that the 1/f part of the metric still satisfies the same 



equations as it does in general relativity. 
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the properties of the stationary, flat-spacetime Klein-Gordon equation in global inertial co- 
ordinates with delta-function source (i.e., the Green's function for the Poisson equation), we 
have A^[<^1 = -levrct"^]; and since c^'^^ is the coefficient of 1/r in equation (125]) . it is in fact just 
the scalar charge q of the body exterior scalar field 0^°^; therefore we take N^'^'^ = — IGvrg. 

Incorporating M and q and rewriting equations f l261l7r|) and ( 131115^ in covarianlo form, 
we have 

E^i^'^ =8n [ MuauMx, z{r))dT (33) 

J 7 



-167r / q6i{x,z{T))dT, (34) 

J 'Y 



as well as 

Mu^Vbu" = q {g"^ + Vfc0(°) (35) 
u'^VaM = -qu^'W a<p'^''\ (36) 



where W a4>^ ' is evaluated on 7. Equations f l5^ and flM]) give the metric and scalar perturba- 
tions produced by the body, showing that they are in fact sourced by the usual point charge 



stress-energy and scalar charge density 



dline 



|8|. Equations (135|) and ( 136|) give the wor 
and mass evolution, and agree with the equations normally given for scalar charges 
Note that the charge q is unconstrained; a separate postulate about the body — such as con- 
stant charge or some other evolution law for q — is required to obtain a deterministic set of 
equations. (Since there is no "conservation law" for g, the body can modify it via internal 
dynamics at will.) Thus the possible motions small bodies in Einstein-scalar theory are 
specified by one free function of time. Equations fl33tl36p give the universal behavior of small 
bodies in Einstein-scalar theory and comprise the results of this subsection. 



B. More General Scalar- Tensor Theories 



The analysis of the preceeding subsection carries over straightforwardly to many more 
general scalar-tensor theories. In fact, the analysis already applies to the majority of scalar- 
tensor theories commonly considered, since these theories have an "Einstein frame" (i.e., a 
field redefinition) in which the matterless Lagrangian (all we ever consider) reduces precisely 

^ Note that spatial components of a tensor Ti in Fermi normal coordinates correspond to projections or- 
thogonal to 7, {5 J' + Uau'')Tf,. 
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to equation (fT5|) . However, suppose that a "Jordan frame" derivation is desired, or that one 
considers a scalar-tensor theory with no Einstein frame. In fact, the analysis still carries 
over to these cases with essentially no modification. Specifically, consider an action 



such that C{g"'^,(j)) is diffomorphism-covariantO and such that E^l = (—g)^^^'^5S/5g"''' and 
^[0] = [—g)^^/'^5S / Scj) are second-order (local) differential operators on {gah^ 0}- The Bianchi 
identity is again given by equation flTIJl) . where now E^^l and EI"^! refer to the new action 
fl37|l . Now adopt the same assumptions for the one-parameter-family {gab{\), 0(A)}, leading 
to equations ( |2TH23|) . Since E^^l and E^'^'^ are assumed second-order, the analysis of appendix 
|A] again gives the effective distributional sources, equations fl26|) and (1271) . Finally, iden- 
tical computations based on the linearized distributional Bianchi identity in Fermi normal 
coordinates give equations ( 13T|) and ( !32i) (as well as N^^ = 0). 

At this point in the treatment of Einstein-scalar theory the parameters Nqq and A^'"'^' 
appearing in equations (!3T|) and (!32l) were interpreted via an analysis of the specific field 
equations for Einstein-scalar theory, where it was found that A^qq = Atccqq = SttM and 
ATM = -levrcl'^] = -167rg, where M and q are the conventional notions of mass and charge. 
For a general theory, the relationship between {A'qq (t), N^'^^{t)} and {c^i/(t, 6, 0), cf'^](t, 6, (f))} 
will depend on the details of the field equations and will in general be more complicated; 
furthermore, there may not be standard notions of mass and charge available. I will simply 
define "mass" M = I/SttNqq and "charge" q = — l/167rA^['^l for a general scalar-tensor 
theory. For any particular theory, one may determine a formula for M and q in terms of c^i, 
and ct*^^ , enabling the calculation of M and q for a particular body from its exterior field via 
the appearance of c^,^ and c^'^] in equations ([HD and (1251) . With M and q incorporated, the 
results for a general a general scalar-tensor theory are ( l33ll36l) . 

Note that one can easily consider theories with multiple scalar fields as well. That is, 
suppose that the Lagrangian of (1371) depends on a whole collection of scalar fields 0/. The 
assumptions are then made for each (pj, and the steps of the derivation proceed apace, with 
copies of equations for each (pj as well as sums over I where appropriate. For example, the 
right-hand-sides of ( 1T9|) . ( 128|) . (l33l) . (l35l) and ( l36l) become sums with one term for each 0/, 
while equations (12111271) and (l3ll) are copied for each 0/. Thus the results are that there is 

i.e., such that £ satisfies ?A*<?!>) = i]uC'{g°'^ ^(t)) for difFemorphisms tp 




(37) 
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a charge associated with each scalar field, and that the force on a body is the sum of the 
ordinary scalar force from each 0/ (and likewise for the mass evolution). Note that if one 
adopts the additional assumption on the field equations for a particular 0/ that g'^ = for all 
bodies (i.e., that bodies do not "produce" this field), then the field 0/ has the interpretation 
of being "non-interacting" (except by gravity) and does not appear in the force law. One 
can add "matter fields" to a theory in this way. 

Of course, it is not at all obvious that our assumptions — argued for in the specific context 
of Einstein-scalar theory (equivalently "Einstein frame" scalar-tensor theory) — will remain 
reasonable for a general theory of the form fl57|) . I now give some examples of more general 
theories in which the assumptions do appear to remain reasonable — that is, in which one still 
expects families of solutions smooth in a and /3 to exist. The first example is the Einstein- 
massive- scalar theory, formed by the addition of a mass term to equation f|T5|l . giving the 
Lagrangiarc 

C = ^[R-2 {g^'Va(pVb(p + r^')] . (38) 

In this case, rather than the Coulomb potential family A/r, the example to keep in mind 
is the Yukawa potential family (A/r)e~^/^. Despite falling off faster than any power of 1/r 
at fixed A, the Yukawa potential family is indeed smooth in a and /3; it is simply I3e~"^^. 
Therefore, the addition of the mass term does not appear to pose any obstacle to smoothness 
in a and (3. Notice that the scaled limit a — )■ gives (3 = 1/f, refiecting appropriate body-like 
falloff in the "buffer zone", even though this falloff does not occur at fixed A. 



A second, more complicated example concerns so-called "chameleon" theories [10|, in 
which non-linear effects, as well as non-minimal coupling to matter, cause a body's exterior 
scalar field — and hence its inferred scalar charge — to depend on 

Although no exact solutions with chameleon behavior are known, approximate solutions that 

have been compared to numerical solutions [l^ show that the exterior field of an isolated 

body is Yukawa if the ambient density is constant (as expected from the linearization of the 

scalar field equation). Since a and /3 near (0,0) refers precisely to the "buffer zone" where 

the density of the external universe would be approximately constant, it seems reasonable to 

expect smoothness to hold here. In fact, this type of argument should work for any theory 

^ Note that the constant £ has dimensions of length (even if G ^ 1); the name "massive" for this theory 

comes from the fact that h/£ would give the mass of excitations of a quantized (f>-fie\d. 
^ A matter field representing ambient density can be included in the action in the manner described in the 



paragraph above that containing p8|) . 
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with a standard kinetic term, since hnearization in the body exterior will give Yukawa. 

The equations of motion for screened bodies in chameleon theory were previously investi- 
gated by Hui, Nicolis, and Stubbs [11], who obtained the non-relativist ic limit of ( 135|) via a 
variant of the original Einstein- Infeld- Hoffman approach. We have rigorously derived the full 
behavior of chameleon bodies, (|35l) and ( !36ll . with no non-relativistic approximations. For 
a Newtonian body with the scalar field coupling usually [lOj considered, the scalar charge 
q corresponds the body's "screened mass". That is, one has a relationship q = q{M,(l)^^^), 
providing deterministic evolution. However, it seems unlikely that any such universal rela- 
tionship will exist for strong-field bodies. The evolution of the scalar chage q would have to 
be determined by other methods in order to take advantage of (!35|) and (!36|) for strong-field 
bodies. 

C. Scaling and Universality 

In this subsection I digress to point out a connection between scaling and universality 
that arises in this work. Returning to the example of Einstein-massive-scalar fl38|) . recall 
that the Yukawa potential example family /3e~"/^ had the scaled limit a— )'0of/3 = l/f, 
which is not the Yukawa potential but the Coulomb potential. The fact that the scaled limit 
gives a field configuration that is not a solution of the theory can be traced to the theory's 
lack of scale invariance (see appendix [B]). In particular, the Lagrangian ( l38i) does not scale 
homogeneously under the rescalings gab ^'^gab and — j- 0. Rather, if one rewrites in terms 
of the barred metric and scalar field, one has 

^ = A'v^ [R - 2 {r'^a^^b^ + X'r'^')] , (39) 

where R is now constructed from cjab- Since the mass term disappears in the A — limit, 
the equations satisfied by g^p and 0*^°-' are the massless Einstein-scalar field equations (a 
fact easily verified at the level of the equations of motion). This explains the appearance of 
the Coulomb potential in the scaled limit of the Yukawa potential example family. 

The fact that the body exterior fields cjp^p and c^^^-* satisfy the massless equations is an 
indication that our results apply only in situations where the mass term can be neglected in 
the buffer zone outside a body. Thus in particular our assumptions require that the body be 
small compared to £. However, this requirement is not suprising or in any sense additional 
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to the basic requirement of the existence of a buffer zone. Because the scale of variation of 
solutions to the massive Klein-Gordin equation is rigidly fixed by I (in that derivatives^ of 
are of order 0/£), the usual requirement that the body be small compared to the scale of 
variation of the external fields in fact implies that the body be small compared with £. This 
implication is conveniently captured in the mathematics of the scaled limit. 

In a general theory (including the higher-rank theories discussed later) the story will be 
analogous: by construction, the scaled limit picks out a scale-invariant "subtheory," which 
must describe the body exterior field approximately for our results to be useful. Thus one 
obtains universal behavior only in the presence of scale-invariance, a situation reminiscent of 
well-known connections between scaling and universality in condensed matter and particle 
physics. Note, however, that the field near the body is not required to be described by a 
scale-invariant theory; our assumptions only refer to f > ^, and furthermore the results 
depend only on f — oo properties of the body exterior fields. Scale-invariance appears 
only in the buffer zone; scale-non-invariant effects (such as "chameleon effects" or those due 
matter fields) are always allowed to act near the body. 

IV. VECTOR-TENSOR THEORIES 

Vector-tensor theories, the most famous of which is Einstein-Maxwell theory, form another 
important class of classical field theories. Here the Lagrangian depends on the metric tensor 
and a vector field A"^, 

S = j d^xC{g^\Aa). (40) 

As in the scalar case we assume that the Lagrangian is diffeomorphism-covariant and the 
field equation operators E^^} = {—gY^'^5S/5g°'^ and £'[^1'' = {—gYf^dS/ 5Aa are second order 
(local) differential operators. A useful example to keep in mind is Einstein-Maxwell theory, 
where (with the normalization of ^\) one has E^^} = Gab-SnT^^^^^ and El^l"^ = -SVaV^'^^]. 
Next we derive the Bianchi identity for a vector theory. Varying f HOj) with repsect to an 
infinitessimal diffeomorphism, we have 

= J rf^xv^|ES(-2V"e') + E^""^^ (fVeA, + A^VaO] , (41) 



This language is slightly sloppy. One really means that scales of variation as measured by the metric, 



such as VTg^^^^V^, are of order 
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for a vector field After integration by parts, the fact that ^'^ is abitrary gives 

As in the scalar case, this identity describes how any "extra" stress-energy E'f^ must be 
non-conserved in the presence of any "extra" charge- current E^^^"" for consistent coupling. 
Note that this identity may also be derived by varying with respect to the upper-index 
A". (However, the appearance will be different when expressed in terms of {—gY^'^SS/Sg"'^ 
computed at fixed A"-, since this quantitity differs from E^^^ (computed at fixed Aa) by terms 
proportional to E^^^"".) We restrict without loss of generality to a lowered-index dynamical 
variable in this section. 

An important special case of vector theories are those whose Lagrangian posesses the 
Maxwell gauge symmetry Aa ^ Aa + Vq'?/'. In this case an anlogous calculation gives 
y^^[A]a = as an identity (describing the requirement that any "extra" charge- current be 
conserved). Thus for theories with the gauge symmetry we have two identities, 

V"eS = i^l^l-^Vi^A] (43) 
Va^I^l'^ = 0. (44) 

The assumptions for a vector field can be motivated by considering the example of 
Einstein-Maxwell theory. Analogously to the scalar case, the sort of behavior we desire 
is the represented by the Coulomb field family ~ A/r. To characterize this type of be- 
havior with a scaled limt, we must define Aa = X~^Aa, so that the scaled limit recovers the 
Coulomb field Aq ~ 1/f. This is also the scaling of Aa that leaves the Einstein-Maxwell 
theory invariant. The analog of (E]) for the components of the vector potential now holds, 
leading again to our assumptions for each component A^ in the coordinates {t, x*). Note that 
our assumptions are on especially strong footing in ordinary electromagnetism, since they 
were in fact shown to hold for the retarded solution of a family of shrinking charge-current 
and stress-energy sources in fiat spacetime [Q]]. I will make these assumptions for a general 
theory of the form (HOl) . 

Therefore the assumptions for this section are the original assumptions of section [Tll with 
the metric components replaced by the pair {g^,y,A^}, which must satisfy E^^^ = and 
jT;[A]a _ g ing^gg^(;[ Qf Eiusteiu's equation. The computation of the small body equations 
of motion proceeds in precise analogy with the scalar case of section lllli That is, define 
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= A^{\ = 0), = dxA^{\)\x=o, and A^^ = \imx-,o A-^{X) (limit at fixed x^^). Then 
the assumptions give series expansions, 

4°)=6[^l(t) + 0(r) (45) 
4^ = 4^^^,^, 0)^ + 0(1) (46) 

A^'^ = bfito) + cjf](to, e, 0)i + O , (47) 

for some b\f^ and c[i*'. The body exterior vector field Ap'^ is seen to be stationary and 
asymptotically flat, confirming its interpretation. Again one finds the effective distributional 
sources at linear order, 

^[.]{i)[^(i),^(i)] = jV^j(t)5(3)(xO (48) 
^(1)] = iv[^]A'(t)5(3)(a;i)^ (49) 

where E^^}^^^ and E'l^K^)'^ are the linearizations of i^jfj and E^^'^^, respectively, off of the 
background {gjj^J , A^/^^}. At this point it makes sense to treat separately those theories 
with the Maxwell gauge symmetry and those without. For those with the symmetry, we 
have the identities (H3l) and (jSj). Employing as usual the linearized, distributional forms 
of these identities in Fermi normal coordinates, we find A^^q' = A^'^' = A"]"^' = as well as 
doN^^^ = SoA-^^^ = and Mai = N^^^^d[oAf. As usual A-^^' and A^^^^ may be interpreted 
by their appearance in the body exterior metric and vector field. In the Einstein-Maxwell 
case discussed above one sees that the usual notions of mass M and charge q are related by 

use this to define q and M for a general theory with 
Maxwell gauge invariance, and rewrite the results covariantly to obtain 

E[f'^ = 8nM I UauM^. z{T))dT (50) 

^[A](l)a _ ^g^^ I u''5i{x,z{T))dT, (5l) 

as well as 

Mv'Vma = qu\2VyaAf), (52) 

where q and M are constants. Thus for theories with the Maxwell gauge symmetry we have 
the usual point particle stress-energy and charge-current, along with the Lorentz force law. 
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For theories without the gauge symmetry, however, the situation is more comphcated. In 
this case we only have the single identity fj42|) . and the Fermi coordinate calculation now 
gives 

= ^ATt^l 4°) (53) 
doNlfJ = -iV[^]^9[,4°) + ^APdoN^^\ (54) 
iViJa. = N^^^^d,,Af + ^5o(ArI^Ur) - ^Af ^^oiVr . (55) 

Since we no longer have A^]q' = A^^^^ = A'^^^"^^ = 0, the parameters q and M no longer suffice 
to characterize the body. In light of equation ( |53l) . which requires A^^'^'j to point along Af^\ 
it seems simplest to introduce a second charge q by N^'^\ = IGnqAf^K Then the results are 
rewritten covariantly as 

= 8^ / {Muau, + 4gP%A)AW) 6,{x, z{T))dT (56) 

J-y 

E[A](l)a ^ ^g^ f fa ^ ^pab^(S>)\ ^^^^^ z{r))dT (57) 

and 

(M - qiA'^A.Y) u'^VaUb = 2{q - gA%,) V[b^]u'^ 

+ P% {2qA'VicAa] + u'Vc (qAdu'^Aa) - A^u''^ ,q] (58) 
u'^VaM = -qA''u\2ViaAb]) + u'^AaU^V^q. (59) 

where P\ = 6\ — projects orthogonally to u", and the superscript (0) on A^a^ has 
been dropped in the last two equations for readability. Thus for theories without the gauge 
symmetry the usual point particle stress-energy and charge-current are not obtained (in 
that the distributional forms are not parallel to u"), the Lorentz force law is not the correct 
force law, and the particle is described by a time-dependent mass M as well and two charges 
q and g, neither of which has an evolution law (just as there was no law for the scalar 
charge in scalar-tensor theory). An example of a commonly studied theory to which these 
equations apply is Einstein-Proca theory. Of course, in most references Proca lagrangian is 
coupled to matter via an interaction term that by itself has the gauge symmetry, so that 
solutions to the Proca-matter system respect charge conservation (but charge conservation 
does not hold as an identity). If one restricts to such matter, presumably one would have 
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q = and u°''Va(l = for bodies made of that matter, whence the usual particle equations of 
electromagnetism f l50ll52|) would be recovered. However, it is not clear that there is reason 
for matter to conserve charge in a vector theory beyond the fact that it is required in the 
most familiar vector theory. It should be emphasized that the Proca theory (and others 
without the gauge symmetry) admit far more general behavior in the motion of bodies than 
does ordinary electromagnetism. This general behavior would give the motion of any matter 
that did not conserve charge, as well as the motion of any non-matter objects (such as black 
holes or "geons" ) that might exist in the theory. 

Note finally that the analysis of this section can be straightforwardly generalized to 
the case of multiple vector fields (or even multiple vector and scalar fields) in the manner 
discussed in section IIIII for scalar-tensor theory. In this case one simply obtains copies of 
equations ( l56l) and ( 157|) for each field, and the right-hand-sides of ( l58l) and ( 159|) are copied 
for each field to form a sum. If any of the fields have Maxwell gauge-invariance, of course, 
the simpler terms from equations ( 13011321) may be used for that field. There may also be 
different gauge symmetries that provide different simplification. An important example is 
non-Abelian gauge theory, where the charges respect "gauge covariant" conservation. More 
precisely, if we label the set of vector fields by and their field equation operators by 
]7;ia ^ (—g)~^/'^5S/5Aj^ (using capital latin indices for "gauge indices") then the gauge 
symmetry of the Lagrangian gives VaE^"' = Yli j f^'^^E^°'A^^ as an identity, where f^'^^ 



are the structure constants as defined in 12| (anti-symmetric on the first two indices). The 
usual Fermi normal coordinate calculation on this identity implies that the "hatted charge" 
vanishes for each body and gives an evolution law for the charges . The small body 
equations of motion are then 

Mu'V.Ua = g'(2V[X])«' (60) 
u'^WaM = J2 f'^'q^KAWu' (61) 

I,J,K 

u^'Vaq^ = J2 f^^Q^^iu". (62) 

The evolution is now fully deterministic on account of the extra symmetry of the Lagrangian. 
Note that if the structure constants are totally anti-symmetric, then M is constant and these 
reduce to "Wong's equations" [3]. 
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V. THE GENERAL CASE 



The procedure used for scalar and vector fields generalizes straightforwardly to higher- 
rank tensor fields. Here we motivate the scaling by analogy with the scalings used before. 
The scalings used for the metric, scalar, and vector fields were all such that the power 
of A cancelled powers of A resulting from the Jacobian of the coordinate transformation to 
scaled coordinates, allowing (|2]) and its scalar and vector analogs to hold. For a general-rank 
tensor field T"'-""^^ j^^, I define the scaled version T"^-''%^ = A"-™r''^-"%^ so that 
the analog of (|2]) holds. Then the usual reasoning leads to the usual assumptions for each 
component of the tensor field, and one proceeds exactly as in the previous sections. I will 
summarize this procedure in the form of a proof of a theorem, below. It is straightforward 
to follow the steps to determine the force law for any theory of particular interest (although 
higher-rank fields are more seldom considered). 

Theorem. Let S be an action in four spacetime dimensions in the sense of appendix E of 
Wald Ji^y such that 1) the Lagrangian C depends differomorphsim-covariantly on the met- 
ric gab and some set of tensor fields {'ipi} (tensor indices suppressed), i.e., (f)*C{gab,'ipi) = 
C{(j)*gab,(f)*'ipi) is satisfied for diffeomorphisms (p; and 2) E^^l = {—g)^^/'^6S / dg""^ and = 
{—g)~^^'^5S/5-ipi are (local) second-order differential operators on {gab,4'i}- Suppose there 
exists a one-parameter-family {(^^^(A), ^/^/(A)} satisfying the analogs of the assumptions of 
section ITWH Then, the worldline 7 (four-velocity u"" and four- acceleration a"") and a cer- 
tain function M defined on 7 satisfy equations of the form Ma'^ = f"" and M°VaM = F , 
where fua = and both f and F are local tensor functions of u"" , a", gab{^)\'y, ''Pi{^)\'y, 
V {ip , and certain tensor fields defined on 7. 

The proof is essentially to follow the steps of the previous sections. Since these steps 
are by now familiar, I will omit some details in the description here (allowing considerable 
savings on notation). Begin by varying S with respect to an infinitessimal diffeomorphism 
to derive the "Bianchi identity" . This results in an expression of the form 

V^i^S = 5^ [ (i? ^ V^,) , + ( Vi? ^ J (63) 



That is, adjoin the components of the "0/ and "0/ to those of the metric and rescaled metric (respectively) 
where they appear, and replace satisfaction of the Einstein equations with satisfaction of the field equations 
E^^l = and = 0. 
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where the notation {A B)^ indicates a sum of terms, each of which consists of the tensor 
product of A and B contracted in some way to yield a dual vector. Note that the explicit 



form for an arbitrary-rank tensor field was worked out in Now derive the effective 

point particle description. By assumption iii), the components of g^^f^^ and tpf'^ are 0(l/r) (r 
near zero). By the assumption of second-order field equations and the analysis of appendix 
lAl the components of e]^^^^^ and E^^^^ as distributions are multiples of 5^(x*); take the 
coefficients to be Njfu{t) and N^{t) (component indices suppressed), respectively. Now apply 
the linearized, distributional form of equation (163|) . By satisfaction of the background field 
equations {E^^^ = E^^^^ = 0), this identity takes the form 

WEl^^W = J2 [{e''~''> V4°))^ + (VE^(I) ) J , (64) 

^ u u 

where here and below V is the derivative operator associated with the background metric. 
We then have 

V^iV^J = Yl [(^' VV^f 1,)^ + (VN^ V^f I,) J . (65) 

^ u u 

Adopting Fermi normal coordinates, one has 

5 (3) (f ) [-^oiVoo + a,N,o + (f ) [N,o - g^] = (66) 
5(3) (f ) [-doNoi + ajNi, + aiNoo - Hi] + di5^'^ (x) [iV,, - Kij] = 0, (67) 

where J-", Qi, Hi, and JCij are local functions of A^^, ipf'''\-y, dipf^l^, ai. Naming A^oo = M, it 
now follows from (166|) and (!67|) that 

doM = GiQi + T (68) 
Mai = doQi - ajJCij + Hi. (69) 

These are the small-body equations of motion for the theory (expressed in Fermi normal 
coordinates). It is straightforward to determine the functions J-", Qi, Hi, and ICij by direct 
calculation for the theory in question (or even the general case; however, the expression 
is not simple). It is also straightforward to relate the parameters M and A^^ to the body 
exterior configuration {gjip , 'ipf'^ } for each particular theory (as done in previous sections), 
enabling their calculation for any particular body. The covariant translation of equations 
(!68|l and (|69|1 proves the theorem. 
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The theorem estabhshes that a simphfied description of motion (in the form of a second- 
order equation for 7]^ is obtained via buffer-zone dynamics in a very large class of theories. 
It is interesting to speculate on the extent to which the class of theories could be enlarged. 
The requirement of second-order field equations seems easiest to relax, since higher-order 
equations would lead only to higher-order delta functions (i.e., derivatives of delta functions) 
appearing in the effective stress-energy, whence our calculations could proceed striaghtfor- 
wardly. However, it is far from obvious that our assumptions on one-parameter-families 
would remain reasonable in the context of higher-order theories, whose solutions may have 
very different propertieso The requirement of a diffeomorphism covariant Lagrangian ap- 
pears difficult if not impossible to relax, since the Bianchi identity plays an essential role 
in determining the motion. At least in our approach, the diffemorphism-symmetry and La- 
grangian formulation are key to obtaining a description of motion from the buffer- zone field 
equations alone. 



VI. SUMMARY 

I have treated the motion of small bodies in classical field theory via the approach of 
{4]. The search for one-parameter-families of solutions representing the exterior field of a 
shrinking body led precisely to the physical assumption of a "buffer zone" — a region far 
enough from the body that its field can be approximated in a multipole series, but close 
enough to the body that the field of the external universe can be approximated in an 
ordinary Taylor series. No assumptions about the body interior are made. In the case of 
second-order metric-based theories following from a diffeomorphism-covariant Lagrangian, I 
derived the force law for scalar and vector fields, and showed that the method works in the 
general-rank case. This provides a rigorous derivation of the small-body force law in many 
classical field theories commonly considered, and shows that that field dynamics outside a 
body determines its motion in a very general class of theories. 
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In some theories the field equations may enforce AI = in which case the equation for 7 may be fower 



than second-order or even triviaL 

Note, however, that many higher-order theories (such as those whose Lagrangian is a function of the Ricci 



scalar) admit second-order formulations, so that the current analysis applies. 
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Appendix A: Delta-function Calculation 

Consider a linear, second-order differential operator L that takes tensors of rank 

n into tensors (L0)^^ of rank m. The adjoint is the linear map from tensors of rank 
m to tensors of rank n defined by jl5| 

for arbitrary if) and (f). The vector s'^ is a multilinear function of 0'»i---<^" VcV'''^'"''") 

and (j)°-'^---°'"^bi...bn _ j^q^ promote L to an operator on distributions in the standard way. 

That is, we define a distributional operator L by 

for smooth test tensors f^^-^-^ of compact support. Now suppose that the metric and 0-field 
have the following expressions in coordinates 

g^.u = V^.u + 0{r) (A3) 
(P^"'-"" = 9, 0) + 0(r°), (A4) 

and further that (j)i^'^---t^" solves the equation at r > 0, i.e., that we have 

(^0)^1-^- =0, r > 0. (A5) 

Here r, 9, (p are related to in the usual way. The distribution [LcpY^---"'^ may be computed 

by 

{Urn = hm / d'x^{L^fr-^"<j>,,,„,^ (A6) 
= hm [ d^x^ {/^— + Vps"} (AT) 
= lim / sin 9d9d(j)dt UpS^, (A8) 

Jr=e 
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where Up is the unit normal to the unit two-sphere in Euchdean three-space, and s'' is 
evaluated on (p and /. The first line follows from flA2l) and the smoothess of ft^^---i^". The 
second line follows from (lAip . The third line follows from (lASp and integration by parts, 
where the boundary term vanishes by the compact support of J^i - a*". The volume element 
on the surface r = e has been replaced with the Minkowski volume element on account of 
the limit e — )■ and the metric form (lA3p . (Note that if different coordinates were chosen 
such that the metric is not Minkowski at r = 0, the volume element would take a different 
form, but it is easy to check that the analysis would still hold.) Now, equation (]A4|) together 
with the properties of imply that s'' may be written 

s'' = iD^^-^-(t, 0,0) (t, r, ^, </)) + 0(l/r) (A9) 
for some Then we have 

{L(l))[f] = Mm [ sin eded(f)dtnpDP'''-''"^ f^, ^Jr=e (AlO) 

= J dtf,,,„,Jr=oJ smeded(l)npDf"''-'- (All) 

^ I dtN''^-''-{t)U„„,Jt,x^ = 0) (A12) 

where the second step follows from the smoothness of f^^-'^rn^ g^^^j g^gp simply defines 

the result of the angular integral to be N^'^-'^"^{t). Since the test function is evaluated at 
= 0, this expression shows that distribution [LcpY'^---^"' is proportional to the spatial delta 
function 5*^'^^(x); i.e., we have the desird result 

{LcPY^-"^ = iV'^^-^'"(t)5(3)(f). (A13) 

Note that an explicit formula for N'^^-'^'^{t) in terms of Cfj^-^,,,p„{t,6,(j)) can be determined 
for any particular differential operator L by following the steps of this computation explic- 
itly. This formula will in general involve angular integrals of C^i...^„ and its first angular 
derivatives. 



Appendix B: Scale Invariance 



In this appendix I provide a definition of scale invariance and compare with other defini- 
tions commonly given. Tensor indices are suppressed throughout this section. It is conve- 
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nient to restrict to theories that follow from a Lagrangian, 

S{g,^j)= I d^xC{g,^i), (Bl) 



although the metric need not be dynamical and there may exist other "background struc- 
ture" . (That is, the Lagrangian is not required to depend diffeomorphism-covariantly on g 
and ipi.) The theory is scale invariant if the action scales homogeneously under a scaling of 
the metric and other fields; more precisely, if exist numbers {Pi} and n such that 

S{\^g,X'''ij) = rS{g,iJi), (B2) 

for numbers A. This property implies that the equations of motion for the rescaled fields are 
identical to the equations of motion for the original fields, so that there is no preferred scale 
for lengths (as measured by g) or field values (0) in the theory. The choice of for the 
scaling of g is conventional, and gives A the interpretation of a length (since the line element 
scales as A^). Note that our definition of ^ = X^^g in the body of the paper is consistent 
with this convention of "(7 — > A^f?". 

An alternative definition of scale-invariance is often given in the context of theories spec- 
ified by partial differential equations in coordinates (without the introduction of a metric). 
In this case one rescales the coordinates — )■ Ax^ and asks if a rescaling of the fields 
can restore the original partial differential equations. Such theories can often be rewritten 
diffeomorphism-covariantly via the introduction of a fiat metric gab- Changing the coordi- 
nates of the original partial differential equation then corresponds to applying the diffeo- 
morphism a associated with coordinate rescaling to all fields except the metric. That is, to 
check if a diffeormpism-covariant equation E[g,il)i\ = is scale-invariant according to the 
non-tensorial definition, one asks if (j^E[g, ipj] = X^'Elg, X^'a^^ipi]. Applying a"^ to both sides 
and using the diffeomorphism-covariance of E, this becomes E[g,iljj] = X'^E[a^lg, X^'tpi]. 
But since g is fiat, (J^\g = X^g, and this reduces to the definition involving rescaling the 
metric. This is why one says that rescaling the metric is a curved-spacetime generalization 
of rescaling the coordinates [l^. 

In the context of quantum field theory a second alternative definition of classical scale- 
invariance is often given by requiring that the action be left invariant under under the 
rescalings, rather than just scale homogeneously This corresponds to our definition with 
the additional demand that n = 0. Since the overall scaling of the action (i.e., the value of n) 
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has no effect on the classical equations of motion, this notion in fact removes as "classically 
scale-invariant" many classical theories with no preferred scale. For example, vacuum general 
relativity is scale-invariant in our sense, having no preferred length scale; however, vacuum 
general relativity is not classically scale invariant according to the definition used in quantum 
field theory. Another example is massless 0" {n ^ 2) in flat spacetime (scale-invariant in 
our sense); only 0^ is classically scale- invariant in the sense used in quantum field theory. 

The difference between the definitions can be further elucidated with reference to the 
well-known fact that a scale-invariant action contains only coupling constants that are di- 
mensionless in particle physics units (c = l, ^=1, Gt^I). Adopting the viewpoint that 
masses and lengths are fundamentally different (but that time intervals and lengths are not), 
it is convenient to work in special relativity units (c = l,h ^ 1,G ^ 1), where the statement 
of dimensionlessness in particle physics units becomes the property of having equal mass 
and length dimension, so that in particular one cannot construct a length using h (which 
has dimensions of mass times length). For example, the constant A of A(j)'^ has dimensions 
of mass times length so that no length can be constructed with h. On the other hand, the 
constant G of general relativity has dimensions of mass over length, so that h can be used to 
construct a length (called the Planck length). While our notion of classical scale-invariance 
implies that no lengths can be constructed from the coupling constants alone, the quantum 
field theory notion of classical scale-invariance places the further restriction that no lengths 
can be constructed even when one is allowed to use h in addition. 
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O ■ Abstract 
(N , 

I show how prior work with R. Wald on geodesic motion in general relativity can be generalized 
to classical field theories of a metric and other tensor fields on four-dimensional spacetime that 1) 
ly-^ ■ are second-order and 2) follow from a diffeomorphism-covariant Lagrangian. The approach is to 

consider a one-parameter-family of solutions to the field equations satisfying certain assumptions 
designed to refiect the existence of a body whose size, mass, and various charges are simultaneously 



, scaled to zero. (That such solutions exist places a further restriction on the class of theories to 

. which our results apply.) Assumptions are made only on the spacetime region outside of the body, 

> 

lO ' so that the results apply independent of the body's composition (and, e.g., black holes are allowed). 

' The worldline "left behind" by the shrinking, disappearing body is interpreted as its lowest-order 

in 

^Sj . motion. An equation for this worldline follows from the "Bianchi identity" for the theory, without 

o . 

, use of any properties of the field equations beyond their being second-order. The form of the 
force law for a theory therefore depends only on the ranks of its various tensor fields; the detailed 



^ ■ properties of the field equations are relevant only for determining the charges for a particular body 

H : 

(which are the "monopoles" of its exterior fields in a suitable limiting sense). I explicitly derive the 
force law (and mass-evolution law) in the case of scalar and vector fields, and give the recipe in the 
higher-rank case. Note that the vector force law is quite complicated, simplifying to the Lorentz 
force law only in the presence of the Maxwell gauge symmetry. Example applications of the results 
are the motion of "chameleon" bodies beyond the Newtonian limit, and the motion of bodies in 
(classical) non-Abelian gauge theory. I also make some comments on the role that scaling plays in 
the appearance of universality in the motion of bodies. 
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I. INTRODUCTION 



In special relativity, a non-interacting body moves in a straight line. Therefore, it is not 
surprising that in general relativity an "infinitesimal test body" (i.e., a body small enough 
that the curvature of the external universe can be neglected, and weakly-gravitating enough 
that curvature it generates can be neglected) will move locally in a straight line, i.e., it 
will follow a geodesic. But from this perspective it does seem quite surprising that strong- 
field bodies like neutron stars and black holes in fact also move on geodesies (in the limit 
of small size). After all, no matter how small or light such a body, the local spacetime 
metric will differ significantly from that of fiat spacetime, and one would therefore expect 
that nonlinear gravitational dynamics — certainly not special relativity — would principally 
determine its motion. Furthermore, since the metrics of different strong-field bodies will 
differ greatly from each other, one would perhaps expect there to be no universal law for 
the motion of strong-field bodies at all. Indeed, the natural assumption would seem to be 
that the motion of a strong-field body depends in detail upon its composition. 

This expectation is incorrect for a very counter-intuitive reason: in general relativity, 
the motion of a small body is in fact completely determined by field dynamics outside of 
the body. This surprising fact was first demonstrated by Einstein, Infeld and Hoffman 
and has become the foundation of a more modern approach to motion termed "matched 
asymptotic expansions" ^ (see also js-G]). The basic physical requirement of this line of 
work is the existence of a region (the "buffer zone") sufficiently far from the body that the 
body field may be approximated as a multipole series, yet sufficiently close to the body 
that the field of the external universe may be approximated in an ordinary Taylor series. 
The vacuum gravitational dynamics taking place in this region then suffice to determine the 
motion. 

A primary purpose of this paper is to determine to what extent this conclusion general- 
izes to other classical field theories. To investigate this question I generalize the approach 
taken in j4| to deriving geodesic motion in general relativity^] In the formalism of 4| a 
small body is characterized by a one-parameter-family of solutions to the vacuum Einstein 
equation describing the region outside of a body that shrinks to zero size and mass with the 
perturbation parameter, A. A family with such behavior is considered by demanding the 



-'^ I do not treat self- force corrections, which were the primary focus of [J]. 
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existence of a second, "scaled" limit wherein the coordinates and metric are rescaled such 
the body is held at fixed size and mass. At A = in the original limit the body disappears, 
leaving behind a smooth spacetime with a preferred worldline, 7, picked out; this worldline 
is interpreted as the lowest-order perturbative motion of the body. We showed that 7 must 
be a geodesic by applying the Bianchi identity to an effective point particle description that 
(remarkably) emerges at first order in A. In this paper I generalize the approach to theories 
that 1) follow from a diffeomorphism-covariant Lagrangian, ensuring a "Bianchi identity" 
and 2) have second-order field equations. 



For the above class of theories the method of j4| gives an equation for 7 that depends 
only on "buffer zone" field properties, showing that the Einstein-Infeld-Hoffman idea remains 
correct in a more general context. More specifically, the equation involves, in addition to the 
value and first-derivative of the external fields at the location of the body, various charges 
(understood to include mass as the charge associated with the metric) that are determined 
from the body's fields in the scaled limit (they are "field monopoles"). The results rely only 
on properties 1) and 2) above and are therefore surprisingly independent of the details of the 
theory. In particular, the force law depends only on the form of the Bianchi identity, which 
in turn depends only on the ranks of the tensor fields considered (although extra identities 
following from gauge symmetry can greatly simplify the results). Therefore, the expression 
for the force in terms of the charges and external field values is in fact identical across 
theories with the same types of tensor fields and the same gauge symmetries. However, the 
charges associated with a particular body composition will differ in different theories, since 
the relationship between a given source and the field monopoles it generates will depend 
on the field equations. In this way different theories will make different predictions for the 
motion of the "same body," even when the force law is identical. 

In interpreting the results it is useful to distinguish varying degrees of "universality" in the 
motion of small bodies. In the case of general relativity, all small bodies move on geodesies, 
so that their internal structure is completely irrelevant to their motion. In Einstein-Maxwell 
theory, a single number characterizing the body (the charge-to-mass ratio) determines how it 
will move, so that the internal structure is minimally relevant. In scalar-tensor theory, a free 
function of time (the charge-to-mass ratio of the non-conserved charge) specifies the motion 
of a body, so that the internal structure is somewhat relevant. In higher-rank theories a 
finite number of free functions of time characterize the motion of a body. Of these results 



only geodesic motion in general relativity is truly universal in that it applies to all bodies; 
however, I will refer to all of the above results as "universal behavior in motion" , since the 
information required to determine the motion of a small body is reduced from the complete 
description of the body to the knowledge of a finite number of parameters at each time. To 
adopt the language of condensed matter physics, there are thus large "universality classes" 
of small bodies that move in the same way. 

The content of this paper is as follows. In section |IT1 1 summarize the formalism of {4] 
to derive geodesic motion in general relativity. In section IIIII I generalize the formalism to 
Einstein-scalar and then more general scalar-tensor theories, deriving the scalar force law. 
Note that mass evolution always occurs, and the scalar charge evolution is unconstrained. I 
discuss the results in the context of specific scalar-tensor theories and comment on scaling 
and universality. In section IIVI I apply the formalism to vector-tensor theories to derive the 
vector force law. This surprisingly complicated equation simplifies to the Lorentz force law 
in theories with the Maxwell gauge symmetry. I also derive the simplified force law in the 
case of non-Abelian gauge theory. Finally in section |V] I give the proof that universality in 
motion is achieved via buffer zone dynamics for tensor fields of arbitrary rank. A definition 
and disambiguation of scale-invariance is given in an appendix. 

I use the conventions of Wald and work in units where G = c = 1. Early- alphabet Latin 
indices a,b,... are abstract spacetime indices, while Greek indices /i, i/, ... give tensor com- 
ponents in a coordinate system. When working in coordinates {t,x^,x'^,x^), mid-alphabet 
Latin indices i,j, ... denote spatial components 1 — 3, while a zero denotes the time com- 
ponent t. Mid-alphabet capital Latin indices I, J, ... label members of a collection of tensor 
fields. 



II. REVIEW OF FORMALISM: MOTION IN GENERAL RELATIVITY 



In this section I review the derivation of geodesic motion given in ^ . While the treatment 



here is self-contained, the reader is referred to j4| for more details and significantly more 
motivation. Note that many of the arguments given here will hold identically or analogously 
for the more general theories treated in later sections, in which case those arguments will 
not be repeated. 

The basic approach to motion is to formalize the notion of a "small body" by considering 



a one-parameter-family of solutions to Einstein's equation that contains a body that shrinks 
to zero size with the parameter A. While no universal behavior in motion (nor even any 
definition of "position"G) can possibly be obtained at any finite A, in the limit A — t- one can 
hope for a simplified description, whose observables will then approximate observables at 
small but finite A fl The task is therefore to develop assumptions on a one-parameter-family 
to the effect that it contains a body shrinking to zero size. The first realization is that the 
body must also shrink to zero mass, since (roughly speaking) no body can be smaller than 
its Schwarzschild radius. The body will thus disappear in the limit, but it will leave behind 
a preferred worldline, 7, characterizing its motion. Our method of considering such a body 
is essentially to demand that if we zoom in on the presumed shrinking body, then a body is 
recovered. This zooming process is accomplished via the notion of a scaled limit, defined as 
follows. Consider a one-parameter-family of metrics Qabi^), whose metric components gfiui^) 
are given in some particular coordinates (t, x*). Introduce both a rescaled metric (jabW = 
^ "^QabW and, for a particular time to, rescaled coordinates (t,x*) = {{t — to)/\,xy\). 
Then, the scaled limit is given by 

gjikto) = \iT^Jtip{^;to), (1) 

A — ^(J 

where the limit is taken at fixed scaled coordinate. In this notation, the bar on the "(yf" 
indicates that the rescaled metric is being considered, while the bars on the coordinate com- 
ponent indices "/x" and "z/" indicate that the components of cjab in the rescaled coordinates 
are being considered. (I will continue to adopt this notation throughout the paperO) 

This limit has the interpretation of "zooming in" because a fixed-x^ observer moves ever 
closer to the shrinking body while the rescaled metric keeps distances finite. A simple 
example to keep in mind is the family of Schwarzschild deSitter metrics of mass parameter A 
^, which clearly contains a shrinking body of the sort we want to consider. In the ordinary 
limit the body disappears, leaving behind the background spacetime of deSitter. But in 
the scaled limit the background "disappears", leaving behind the Schwarzschild metric for 



all to- For all one-parameter-families in this paper, I will refer to gjjj = g^v{\ = 0) 



as 



^ Consider, for example, the impossible task of assigning a center of mass position to a black hole. 

^ Of course, the value of A at which the physical spacetime is embedded into the one-parameter-family is 

arbitrary. What matters for the application of the simplified description is that corrections to the relevant 

observables are numerically small. 
^ Note also that coordinate indices will always refer to the original cartesian-like (t, x'), even if the coordinate 

components are being viewed as functions of other variables, such as spherical coordinates. 
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the "background metric" and to (jp'pito) as the "body exterior metric" (no interior is ever 
considered). The existence of both original and scaled limits reflects the presence of a body 
with the appropriate scaling, and will be our assumptions i) and ii), below. 

An additional assumption is required. To arrive at this assumption, note that the rescaling 
of the metric by effectively cancels powers of that arise in changing to the rescaled 
coordinates, so that one has the simple formula 

^/ip(A; to; t, x') = g^y{\] t = to + At, x' = Ax'), (2) 

relating barred components of the barred metric, g^p, to corresponding unbarred coordinates 
of the unbarred metric, g^^u- That is, one simply "plugs in" t = to + At and = Ax* to 
compute cjpjy from g^^. This formula shows that if we regard g^u as a function of new variables 
a = r,l3 = X/r,6,(j) (at fixed t, and with r, 6, defined relative to x* in the usual way), then 
the scaled limit ([1]) is given by the limit a — ?■ at fixed /3 of the original metric components 
gfj_u. Similarly, the original limit is given by the limit /3 — )■ at fixed a. Demanding that 
both limits exist is thus the statement of separate continuity in a and /3. A natural extension 
(argued for at length in |j] on the grounds that it excludes certain pathological behavior) 
is to demand joint continuity and in fact joint smoothness in a and /3 (although only 
is required here). This will be our assumption iii), below. Note that the electromagnetic 
analog of iii) has been shown to hold for a family of shrinking charge-current sources in flat 
spacetime 8|. 

The assumptions of this section for the metric family gab{X) are that there exists coordi- 
nates (t,x*) at each A > such that the following hold. 

• i) For r > A-R for some constant R, the metric components g^v{X) satisfy the vacuum 
Einstein equation and are smooth functions of (t,x*. A). The worldline, 7, defined by 
A = X* = is timelike. 

• ii) The scaled metric components gp,u(to; A) are smooth functions of (A, t, x*) for f > R. 

• iii) The metric components g^i, are smooth functions of {a, (3) at (0, 0) for fixed (t, 6', 0). 

Assumption i) establishes our domain r > XR and provides the requisite smoothness for 
perturbation theory on that domain. It also lays the groundwork for the interpretation of 
the domain as the exterior of a shrinking body by taking 7 to be timelike. Assumption 



ii) establishes this interpretation according to the ideas of the scaled limit, and assumption 

iii) adds additional "uniformity" properties j^j. Note, however, that this latter assumption 
has an important physical consequence (i.e., it places an important restriction on the type 
of spacetime for which our approximate results will be useful). It requires that there be a 
spatial region both far enough from the body that its field can be approximated in a series 
in inverse powers of distance (corrections in /3 near zero), and close enough to the body 
that the field of the external universe can be approximated as a series in positive powers of 
distance (corrections in a near zero). Therefore by seeking one-parameter-families containing 
a shrinking body, we in fact end up with a mathematically precise version of the usual "buffer 
zone" assumption of the Einstein-Infeld-Hoffman approach and its descendants. No further 
assumptions beyond i), ii), iii) (and their analogs for other theories) are made in this paper. 

Smoothness in a and /3 allows us to Taylor expand in these variables to any finite order. 
However, to derive geodesic motion we in fact require only a single derivative in /3, 

g^,{\ t, r, 0) = 6^,(t, 0) + c^,(t, 0)/3 + 0{a) + 0{f) 

= b^,{t, e, <P) + c,,{t, e,<P)^ + 0{r) + O 0^ (3) 

where 0{r) near zero is at fixed A/r, and 0{X/r) near zero is at fixed r. Sorting into powers 
of A and r, we have 

g,,{\, t, r, 6, 0) = b,,{t, 9, 0) + 0(r) + A (c^,{t, ^, 0)^ + 0{r')\ + 0{\^), r > (4) 



where the order symbols are for small r and A. From this expression it is easy to read off 
series expressions for the background metric gj^J and linear perturbations g^^J = dxg^u\x=o, 

g^^Xt,x') = b,At) + 0{r), r>0 (5) 

g^Kt, x') = c,,{t, 9,<P)^ + 0(1), r > 0. (6) 



Since gj^) is assumed smooth everywhere, cannot depend on angles and we have written 
h^yit). (This quantity is usually taken to be rj^y by coordinate choice.) We are also interested 
in the consequences of equation ([3]) in the scaled limit. Using equation ([2]), we have 

gj,^{X- to; i, f, 9, 0) = b^,{to + At, 9, 0) + c^,(to + At, 9,(f))^ + 0(Af) + O (J^^ , (7) 

so that the limit A — > at fixed (t, x*) gives 

g^kto; x') = b^,{to) + c^uito, ^,0)^ + , (8) 
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and the body exterior metric g^J is seen to be stationary (independent of i) and asymptot- 
ically flat (constant as f — > oo). (Stationarity follows from smoothness of g^uW in ^0 This 
supports the idea that cjp'p characterizes the body exterior as it would appear in isolation at 
time to- Equations (E]), ([6]), and ([H]) are key consequences of our assumptions. 

Geodesic motion is now derived as follows. By assumption i), g^y{\) is smooth in A at 
A = for r > 0. Thus (yf^l] satisfies the vacuum linearized Einstein equation about g^^\ 

G'iim = ^, r>0. (9) 

Now regard g\iv &s a distribution on defined on the background spacetime including at 
r = 0, which is possible because its "most singular" behavior is only 1/r. Since IS a 
second-order, linear partial differential operator, it follows from equations ([6]) and ([9]) via 
the analysis of appendix |X] that, distributionally, we have 

GfM'^] = N,,{t)5^'\x^), (10) 

for some N^j^y defined on the worldline = 0. (This result is analogous to the well-known 
fact that V^(l/r) = —Att6^'^\x^). If the explicit form of Gji] is used, a formula for may be 
obtained for N^y in terms of angle averages of c^,^ and its first angular derivatives.) Thus, 
an effective distributional stress-energy of l/87rA^^,y5^(x*) has emerged at first-order in per- 
turbation theory, supported on the worldline 7. This is remarkable, given that any true 
stress-energy associated with the body is confined to r < A-R and excluded from considera- 
tion; and furthermore, the body need not be "made" of stress-energy at all (as in the case 
of a black hole). 

The strategy is now to apply "conservation" to the "stress-energy" . That is, because the 
linearized Bianchi identity '^"'G^^^[g^^\ = holds as an identity on all sufficiently smooth 
9ab (^*^^ necessarily satisfying the linearized Einstein equation), the distributional linearized 
Bianchi identity also holds as an identity on distributional g'^^^ , and we must have 

W {N^y{t)S^'\x')) = (11) 

in the distributional sense. Here Va is the derivative operator associated with the background 
metric g^^J . The consequences of this equation can be determined in a variety of ways. I 
will proceed by adopting the specific coordinate choice of Fermi normal coordinates (see, 
e.g., {3I) for the background metric gj^J . On the worldline = 0, the metric components 
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are Minkowski (g^^ = r]^y), and the ChristofFel symbols are given by Fqq = Fq^ = Oj, where 
tti are the spatial components of the four-acceleration of the worldline (the time component 
Oo is zero). We then have 

(iV^o^^'H^)) = <^^'H^) [-^oiVoo + a.Ar,o] + d,5^''\x) [N,,] (12) 
{N^i5^^^ (x)) = 5(3) (f ) [-doNoi + a.Nij + a.iVoo] + di5^^^ (f ) [N^j] , (13) 

with repeated spatial indices summed. The coefficients of 5^'^\x) and di5'^'^\x) must sep- 
arately vanish, giving A^jo = = 0, as well as doNoQ = and ajA^oo = 0, i.e., geodesic 
motion when A^oo is non-zero. 

We can interpret Nqq through its appearance in equation (fTOj) . Since it multiplies the delta 
function, it is clear that A^oo will determine the singular behavior of the metric perturbation, 
i.e., it will determine the coefficient c^,y in equation Furthermore, if one expands the 
background metric according to (jS]), it is clear that only the constant term b^i, (which here 
equals 77^,^ by coordinate choice) is relevant for the determination of c^,y via ffTOj) . Thus, 
we may compute this coefficient by using the stationary linearized Einstein equation off of 
flat spacetime in global inertial coordinates. The solutions are well known and one obtains 
for the time-time component (which is all we need) that Atccqq = Nqq (so that cqo is in fact 
independent of angles in these coordinates). Observing the appearance of c^u in the series 
for the body exterior metric (jS]), we conclude that A"oo is Svr times the ADM mass of the 
body exterior metric. Therefore we define M = (l/87r)A^oo and refer to this quantity as 
the mass of the body. This explains the role of the requirement that A"oo 7^ for geodesic 
motion to hold: there must actually be a body present in the one-parameter-family for the 
curve to be necessarily geodesic. 

Equation fllUp for the effective stress-energy may be clarified by introducing M and by 
rewriting in covariant form. Since = (1,0) in Fermi normal coordinates, we have from 
Nij = Nio = that Nab = MuaUb- The spatial coordinate delta function becomes a worldline 
integral of the "invariant" four-dimensional delta function 54{x,x') = Thus we 

have 

G'-ablg^'^] = 8vrM [ UauM^, 4r))dT, (14) 
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where the mass M is constant. This equation shows that the metric perturbations for our 
family are in fact sourced by the usual "point particle" stress-energy (see, e.g., 3[). Thus, 
despite the fact that point particles do not make sense in general relativity (9|, we have 



shown that they do emerge as part of a (mathematically rigorous) approximate description 
of the metric of an arbitrary small body. Furthermore, the "particle mass" M is indeed the 
ADM mass of the body (as measured in the scaled limit). 

The results of this section (i.e., the results of sec. IV of j4|) may be summarized as 
follows. Consider a one-parameter-family of spacetimes containing a body whose size and 
mass decrease to zero, according to the stated assumptions. Then, the ADM mass M of the 
body exterior metric is a constant independent of time to? and, if M 7^ 0, the worldline 7 left 
behind after the body disappears is a geodesic of the spacetime g'^^^ left behind. Furthermore, 
the far-field effective description in linearized gravity is that of a point particle of mass M. 
These results show, in essence, that small bodies move on geodesies while keeping their ADM 
mass constant and sourcing linear perturbations reflecting a point particle of that mass. 



III. SCALAR-TENSOR THEORIES 



A simple generalization of general relativity is the addition of a scalar field. I will first 
consider the ordinary Einstein-scalar theory in detail. I will then discuss the general case, 
which in fact follows from the computations already done. Finally I make some comments on 
scaling and universality. In this and later sections it will be convenient to use a Lagrangian 
formulation. I will use the definitions and conventions of appendix E of Wald [zj, except 
that I will denote his fixed volume element e by "d^a;" . 



A. Einstein-scalar theory 

The action for general relativity plus a minimally-coupled massless scalar field is given 

by 

S = j rf^sv^ {R - 2g''^Va(t)^b(f) , (15) 

where R is the Ricci scalar constructed from g„}). I have chosen the relative normalization 
so that the theory reduces to that of Quinn jlo| in the appropriate limit. It is helpful to 
define E^^} = {-g)-^/^6S/6g^'' and E^^^ = {-g)-^/^S/6(j), which evaluate to 

=Gab-2 (Va<pVb<P - Igabg^'VcCpVdci^ (16) 

= Ag'"'VaVb(j). (17) 
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The equations of motion for the Einstein-scalar theory are then simply e]^^ = and E^'^'^ = 0. 

The formalism requires a "Bianchi identity" for this theory. To derive such an identity, 
consider the variation of the action f|T5l) with respect to an infinitessimal diffeomorphism. 
Since the action is diffeomorphism- invariant, the variation must vanish, and one has 

= I d'x^g {ES(-2V"e') + E^^\'^Va<p] , (18) 

where is an arbitrary vector field. Integrating the first term by parts, we derive 

V^i^S = -^i^'^'Vfcc/). (19) 

The field equations were not used in deriving this equation, which therefore holds as an iden- 
tity on all sufficiently smooth {gab, 0}- (This can also be easily checked by direct calculation 
using equations (fT6!) and (fT7|) .) To interpret this identity, note that nonzero values of e]^^ 
and E'^ would normally be interpreted as stress-energy and scalar charge density (respec- 
tively) associated with some matter field. This equation gives the precise non-conservation 
the matter stress-energy in terms of the matter scalar charge density necessary for consis- 
tent coupling of that matter to the Einstein-scalar theory. Although we will always impose 
^[3] = and = at finite A, non-zero values will emerge in the linearized, distribu- 
tional description (analogously to equation ( |T0|) in general relativity), reflecting an effective 
stress-energy and scalar charge of the body. 

We now seek to generalize the assumptions used in general relativity to the Einstein- 
scalar theory. The main requirement is to take the scalar charge to zero along with the size 
and the mass, in order to keep the energy in the field finite. Thus we seek a one-paramater 
family with scalar field behavior like ~ A/r. To characterize this by the existence of a 
scaled limit, the appropriate rescaling (after changing to scaled coordinates) is simply (p = (p. 
(That is, no rescaling is required; however, we still define (f) for notational consistency.) Then 
for an arbitrary family we define the scaled limit as in ([T]), 

(^(°) = lim(^(A), (20) 

A— >0 

where the limit is taken at fixed scaled coordinate. I will refer to 0'-°^ as the body exterior 
scalar field in analogy with the body exterior metric cjp^p . 

One can now follow the same path of reasoning as in section [Til leading one to assume 
the existence of original (i) and scaled (ii) limits, as well as the uniformity condition (iii). In 
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other words, the appropriate assumptions for this section are those of HIl where the metric 
gabW is replaced by the pair {(^^^(A), 0(A)}, and the required equations are not Einstein's 
equation but the Einstein-scalar equations eI^^ = and E^"^^ = 0. (The same coordinates 
(t, X*) and hence worldline 7 are used for the metric and scalar field.) Note that the new 
assumption iii) will imply the existence of a "buffer zone" for the scalar field as well as the 
metric. 

The steps of the derivation of motion now follow completely analogously. The analog of 
(|3]) holds for the scalar field, which leads to the analogs of ([5]), ([6]) and ([8]), given by 

0(0) =6M(t) + 0(r) (21) 
0(1) =c[<^](t,^,0)l + O(l) (22) 



0(°) = b^^\to) + c[^l(to, e,<p)- + (^) (23) 



for some fo''^' and c^'^l (Here (j)^^'^ = 9a0|a=o is the scalar field perturbation, so that one has 
0(A) = 0(0) + A0(i) + 0(A2) for r > 0.) The body exterior scalar field 0(o) is seen to be 
stationary and to have smooth behavior in 1/f as f — )■ 00, supporting the name we have 
given it. Note, however, that in contrast to the case of the metric, the body exterior scalar 
field retains some "memory" of the external universe, since fe'"^' contains physical information 
about 0(0). As before, the perturbations 0(i) and g'^J satisfy the linearized field equations 
at r > 0, 

E^P[g^'\<j^^'^] = 0, r>0 (24) 

= 0, r>0 (25) 

where -E^^t.'^^^ denote the linearizations of E^f^ and E^'^\ respectively, off of the 

background {g^^\ 0(o)}. Since the field equations are second order in gab and in 0, the linear 
operators e'^^^^'^ and E^^^^-^'^ are also second order in g'^J and (p^-^K Thus, as in the case of pure 
gravity, the "1/r behavior" of equations ([6]) and fl22l) implies via the analysis of appendix El 
that, distributionally, we have 

El!n^'\<l>^'^] = Nl^}{t)6^'\x^) (26) 
^M(i)[^(i)^0(i)] = ivM(t)5(3)(a;i)^ (27) 

for some Njfj and N^'^h As discussed above, a non-zero value of E^^u would normally be 
interpreted as a matter stress-energy source, whereas a non-zero value of E^'^'^ would be 
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interpreted as a scalar charge density source. Thus, despite the fact that no matter stress- 
energy or scalar charge density was considered, distributional descriptions of both have arisen 
effectively in perturbation theory. We can take advantage of this remarkable occurrence, as 
before, by using the "Bianchi identity" for the theory. That is, since the linearization of 
equation V^jfj^^^ = -1/2E['^](i) V^0(°), holds as an identity on all sufficiently smooth 
{gf^v, 0^^)}, it must hold as an identity on distributions as well. Thus we must have 



(28) 



in the distributional sense. Here Va is the derivative operator associated with the background 
metric g^^J . Adopting Fermi normal coordinates as before, this becomes 








5(3) (f) 
5(3) (f) 



N 



[9] 



iO 



+ 9^5(3) (f) 



[9] 



from which we determine 



[g] 

ifi 



as well as 



(29) 

(30) 

(31) 
(32) 



We can interpret A^^qq (t) and A^['^](t) as follows. The formula for E'fjj, equation ( |T6i) . demon- 



Is] 



strates that the linearization, eI^^^^^ will depend on second-derivatives of g^^l' only through 
the linearized Einstein tensor G^^ij , and will contain no second-derivatives of 0^^) . Therefore, 
the identical arguments from th^ case of pure gravity carry over, and we have (l/47r)Ai"QQ = 



(1) 



Coo = 2M, where M is the masj^l of the body exterior metric at time to = t. Similarly, from 
the properties of the stationary, flat-spacetime Klein-Gordon equation in global inertial co- 
ordinates with delta-function source (i.e., the Green's function for the Poisson equation), we 
have A^l"^! = -167rc[<^]; and since ct'^^ is the coefficient of 1/r in equation (123|) . it is in fact just 
the scalar charge q of the body exterior scalar field 0^°); therefore we take N^'^'^ = —IGirq. 
Incorporating M and q and rewriting equations ( 126H27|) and ( 13TH32|) in covarianlo form. 



^ Applying the usual notion of mass to Einstein-scalar theory makes sense because the scalar field stress- 
energy is quadratic in first-derivatives of (f), so that the 1/f part of the metric still satisfies the same 

equations as it does in general relativity. 
^ Note that spatial components of a tensor Ti in Fermi normal coordinates correspond to projections or- 



thogonal to 7, {Sj' + Uau'')Tf,. 
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we have 

Ei^J(^) = StT /" MUaUb5i{x,z{T))dT (33) 

^mi) = _i5^ f qS^(x, ziT))dT, (34) 

as well as 

Mu^Vbu'' = q {g"^ + u^u^) Vfe0(°) (35) 

w'^V.M = -gn^'V^^^o), (36) 
where Va0*-°^ is evaluated on 7. Equations f p3|) and flMj) give the metric and scalar perturba- 



tions produced by the body, showing that t 
stress-energy and scalar charge density jsf 



ley are in fact sourced by the usual point charge 



10[. Equations fl5^ and f l5B]) give the worldline 
and mass evolution, and agree with the equations normally given for scalar charges js], lo| . 
Note that the charge q is unconstrained; a separate postulate about the body — such as con- 
stant charge or some other evolution law for q — is required to obtain a deterministic set of 
equations. (Since there is no "conservation law" for g, the body can modify it via internal 
dynamics at will.) Thus the possible motions small bodies in Einstein-scalar theory are 
specified by one free function of time. Equations fl33ti36p give the universal behavior of small 
bodies in Einstein-scalar theory and comprise the results of this subsection. 



B. More General Scalar- Tensor Theories 



The analysis of the preceding subsection carries over straightforwardly to many more 
general scalar-tensor theories. In fact, the analysis already applies to the majority of scalar- 
tensor theories commonly considered, since these theories have an "Einstein frame" (i.e., a 
field redefinition) in which the matterless Lagrangian (all we ever consider) reduces precisely 
to equation f|T5|) . However, suppose that a "Jordan frame" derivation is desired, or that one 
considers a scalar-tensor theory with no Einstein frame. In fact, the analysis still carries 
over to these cases with essentially no modification. Specifically, consider an action 

S = j d^xC{g''\ 0) (37) 

such that C{g"'^, (p) is diffeomorphism-covarian and such that e[^^ = {-g)-^/^5S / dg"^ and 
j^[4>] = (^—g'j^^/'^§S/5(j) are second-order (local) differential operators on {gab, <!>}■ The Bianchi 
i.e., such that C satisfies C{4'*g°''^ ,i'*4>) = '4'*^{9°'^ :4>) for diffemorphisms '0 
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identity is again given by equation f lT9|) . where now ^2 and El'^l refer to the new action 
( 137|) . Now adopt the same assumptions for the one-parameter-family {gabW, (PW}: leading 
to equations f l2TH23|) . Since e]^^ and E^^^ are assumed second-order, the analysis of appendix 
R] again gives the effective distributional sources, equations (l26ll and ( 1271) . Finally, iden- 
tical computations based on the linearized distributional Bianchi identity in Fermi normal 
coordinates give equations (EI]) and ([32]) (as well as A^^^^ = 0). 

At this point in the treatment of Einstein-scalar theory the parameters A^qq and N^'^^ 
appearing in equations (IHT]) and ([22D were interpreted via an analysis of the specific field 
equations for Einstein-scalar theory, where it was found that A^qq = 471000 = SvrM and 
ATM = -IQnc^'f'] = -IQnq, where M and q are the conventional notions of mass and charge. 
For a general theory, the relationship between {A'qq (t), N^'^^{t)} and {c^!/(t, 0, 0), c['^](t, 6, 0)} 
will depend on the details of the field equations and will in general be more complicated; 
furthermore, there may not be standard notions of mass and charge available. I will simply 
define "mass" M = 1/8ttNqq and "charge" q = — l/167rA^['^l for a general scalar-tensor 
theory. For any particular theory, one may determine a formula for M and q in terms of c^i/ 
and c^'^l , enabling the calculation of M and q for a particular body from its exterior field via 
the appearance of c^^ and c'^'^'^ in equations ([8]) and (l23l) . With M and q incorporated, the 
results for a general scalar-tensor theory are (I33tl36p . 

Note that one can easily consider theories with multiple scalar fields as well. That is, 
suppose that the Lagrangian of ( 1371) depends on a whole collection of scalar fields 0/. The 
assumptions are then made for each 0/, and the steps of the derivation proceed apace, with 
copies of equations for each 0/ as well as sums over I where appropriate. For example, the 
right-hand-sides of (|T9]) . (128]) . (I33l) . (135]) and (136]) become sums with one term for each 0/, 
while equations (12111271) and ( IM]) are copied for each 0/. Thus the results are that there is 
a charge q^ associated with each scalar field, and that the force on a body is the sum of the 
ordinary scalar force from each 0/ (and likewise for the mass evolution). Note that if one 
adopts the additional assumption on the field equations for a particular 0/ that = for all 
bodies (i.e., that bodies do not "produce" this field), then the field (pj has the interpretation 
of being "non-interacting" (except by gravity) and does not appear in the force law. One 
can add "matter fields" to a theory in this way. 

Of course, it is not at all obvious that our assumptions — argued for in the specific context 
of Einstein-scalar theory (equivalently "Einstein frame" scalar-tensor theory) — will remain 
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reasonable for a general theory of the form (1371) . I now give some examples of more general 
theories in which the assumptions do appear to remain reasonable — that is, in which one still 
expects families of solutions smooth in a and /3 to exist. The first example is the Einstein- 
massive- scalar theory, formed by the addition of a mass term to equation flT^ . giving the 
Lagrangiarc 

£ = [i? - 2 ((^'^'Va^Vfe^ + r202)] . (38) 

In this case, rather than the Coulomb potential family A/r, the example to keep in mind 
is the Yukawa potential family (A/r)e~^/^. Despite falling off faster than any power of 1/r 
at fixed A, the Yukawa potential family is indeed smooth in a and /3; it is simply /3e~"/^. 
Therefore, the addition of the mass term does not appear to pose any obstacle to smoothness 
in a and /3. Notice that the scaled limit a — gives /3 = 1/f, reflecting appropriate body-like 
falloff in the "buffer zone", even though this falloff does not occur at fixed A. 

A second, more complicated example concerns so-called "chameleon" theories HI], in 
which non-linear effects, as well as non-minimal coupling to matter, cause a body's exterior 
scalar field — and hence its inferred scalar charge — to depend on 

Although no exact solutions with chameleon behavior are known, approximate solutions 
that have been compared to numerical solutions [ll| show that the exterior field of an 
isolated body becomes Yukawa if the ambient density is constant (as expected from the 
linearization of the scalar field equation). Since a and /3 near (0,0) refers precisely to the 
"buffer zone" where the density of the external universe would be approximately constant, 
it seems reasonable to expect smoothness to hold here. In fact, this type of argument should 
work for any theory with a standard kinetic term, since linearization in the body exterior 
will give Yukawa. 

The equations of motion for screened bodies in chameleon theory were previously investi- 
gated by Hui, Nicolis, and Stubbs [12], who obtained the non-relativistic limit of ( l35l) via a 
variant of the original Einstein- Infeld- Hoffman approach. We have rigorously derived the full 
behavior of chameleon bodies, fl35|) and fl36|) . with no non-relativistic approximations. For 
a Newtonian body with the scalar field coupling usually [11] considered, the scalar charge 
q corresponds the body's "screened mass". That is, one has a relationship q = g(M, 0*^'')), 
^ Note that the constant £ has dimensions of length (even if G ^ 1); the name "massive" for this theory 

comes from the fact that h/£ would give the mass of excitations of a quantized (f>-fie\d. 
^ A matter field representing ambient density can be included in the action in the manner described in the 



paragraph above that containing p8|) . 
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providing deterministic evolution. However, it seems unlikely that any such universal rela- 
tionship will exist for strong-field bodies. The evolution of the scalar charge q would have to 
be determined by other methods in order to take advantage of and (15^ for strong-field 
bodies. 



C. Scaling and Universality 

In this subsection I digress to point out a connection between scaling and universality 
that arises in this work. Returning to the example of Einstein- massive-scalar fl55]) . recall 
that the Yukawa potential example family I3e~°'/^ had the scaled limit a— !'0of/3 = l/f, 
which is not the Yukawa potential but the Coulomb potential. The fact that the scaled limit 
gives a field configuration that is not a solution of the theory can be traced to the theory's 
lack of scale invariance (see appendix [B]) . In particular, the Lagrangian ( 138|) does not scale 
homogeneously under the rescalings gab ^^gab and — )■ 0. Rather, if one rewrites in terms 
of the barred metric and scalar field, one has 

C = A^v^ [R - 2 (^'^'V,0V,0 + , (39) 

where R is now constructed from (jab- Since the mass term disappears in the A — )■ limit, 
the equations satisfied by cjp'p and f^^^-* are the massless Einstein-scalar field equations (a 
fact easily verified at the level of the equations of motion). This explains the appearance of 
the Coulomb potential in the scaled limit of the Yukawa potential example family. 

The fact that the body exterior fields g'^^} and (p'^^^ satisfy the massless equations is an 
indication that our results apply only in situations where the mass term can be neglected in 
the buffer zone outside a body. Thus in particular our assumptions require that the body be 
small compared to I. However, this requirement is not surprising or in any sense additional 
to the basic requirement of the existence of a buffer zone. Because the scale of variation of 
solutions to the massive Klein-Gordin equation is rigidly fixed by £ (in that derivatives^ of 
are of order 0/£), the usual requirement that the body be small compared to the scale of 
variation of the external fields in fact implies that the body be small compared with ^. This 
implication is conveniently captured in the mathematics of the scaled limit. 



This language is slightly sloppy. One really means that scales of variation as measured by the metric, 



such as yJ\g°-^\/a<l)^b(t)\, are of order 
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In a general theory (including the higher-rank theories discussed later) the story will be 
analogous: by construction, the scaled limit picks out a scale-invariant "subtheory," which 
must describe the body exterior field approximately for our results to be useful. Thus one 
obtains universal behavior only in the presence of scale- invariance, a situation reminiscent of 
well-known connections between scaling and universality in condensed matter and particle 
physics. Note, however, that the field near the body is not required to be described by a 
scale-invariant theory; our assumptions only refer to f > ^, and furthermore the results 
depend only on f — )■ oo properties of the body exterior fields. Scale-invariance appears 
only in the buffer zone; scale-non-invariant effects (such as "chameleon effects" or those due 
matter fields) are always allowed to act near the body. 



IV. VECTOR-TENSOR THEORIES 



Vector-tensor theories, the most famous of which is Einstein-Maxwell theory, form another 
important class of classical field theories. Here the Lagrangian depends on the metric tensor 
and a vector field A"-, 

S = j d''xC{g''\Aa). (40) 

As in the scalar case we assume that the Lagrangian is diffeomorphism-covariant and the 
field equation operators E^^^ = {—gY^'^dS / dg""^ and £'1^1'^ = [—gY/'^6S/8Aa are second order 
(local) differential operators. A useful example to keep in mind is Einstein-Maxwell theory, 
where (with the normalization of Q) one has = G^fc - SttT^^*^ and E^^]'' = -SVaV^^'A^l 
Next we derive the Bianchi identity for a vector theory. Varying (140!) with respect to an 
infinitessimal diffeomorphism, we have 

= j d'x^ {ES(-2V"e') + i^'^^" {C^cAa + A^VaO] , (41) 
for a vector field After integration by parts, the fact that is arbitrary gives 

ya^g ^ ^[^]-V[,^] + ^V.^I^'Mfe. (42) 

As in the scalar case, this identity describes how any "extra" stress-energy E'^^' must be 
non-conserved in the presence of any "extra" charge-current E^^^"" for consistent coupling. 
Note that this identity may also be derived by varying with respect to the upper-index 
A"". (However, the appearance will be different when expressed in terms of {—gY^'^SS/Sg"'^ 
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computed at fixed A"", since this quantity differs from e]^^ (computed at fixed Aa) by terms 
proportional to E^^^'^.) We restrict without loss of generality to a lowered-index dynamical 
variable in this section. 

An important special case of vector theories are those whose Lagrangian possesses the 
Maxwell gauge symmetry Aa ^ Aa + Va'ip- In this case an analogous calculation gives 
y^^[A]a = as an identity (describing the requirement that any "extra" charge- current be 
conserved). Thus for theories with the gauge symmetry we have two identities, 

Va^l ^ E^A]a^[aA,] (43) 

Va^l^]'^ = 0. (44) 

The assumptions for a vector field can be motivated by considering the example of 
Einstein-Maxwell theory. Analogously to the scalar case, the sort of behavior we desire 
is represented by the Coulomb field family ~ A/r. To characterize this type of behav- 
ior with a scaled limit, we must define Aa = X^^Aa, so that the scaled limit recovers the 
Coulomb field Aq ~ 1/f. This is also the scaling of Aa that leaves the Einstein-Maxwell 
theory invariant. The analog of (I2j) for the components of the vector potential now holds, 
leading again to our assumptions for each component A^ in the coordinates (t, x*). Note that 
our assumptions are on especially strong footing in ordinary electromagnetism, since they 
were in fact shown to hold for the retarded solution of a family of shrinking charge-current 
and stress-energy sources in fiat spacetime 8!]. I will make these assumptions for a general 



theory of the form (1401) . 

Therefore the assumptions for this section are the original assumptions of section [Tll with 
the metric components replaced by the pair {(7^^,^!^}, which must satisfy e]^^ = and 
j7;[A]a _ g j^stead of Einstein's equation. The computation of the small body equations 
of motion proceeds in precise analogy with the scalar case of section lllli That is, define 
= A^{\ = 0), 4^ = dxA^{X)\x=o, and = limA-,o4(^) O^i^ at fixed x^). Then 
the assumptions give series expansions, 

4)=6jf](t) + 0(r) (45) 
4^ =c[f'(i,^, 0)^ + 0(1) (46) 
4°^ = bfito) + cfito, e, 0)i + 0(^^y (47) 
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for some and c|^^'. The body exterior vector field A^^^ is seen to be stationary and to 
approacfi a constant value as f — oo, confirming its interpretation. Again one finds the 
effective distributional sources at linear order, 

^[.](i)[^(i),^(i)] = iVM(t)5(3)(xO (48) 

where E^^l^^"* and E'l^K^)'^ are the linearizations of i^jfj and E^^^^, respectively, off of the 
background {gj^J , A^/^^}. At this point it makes sense to treat separately those theories 
with the Maxwell gauge symmetry and those without. For those with the symmetry, we 
have the identities (H3l) and (jS]). Employing as usual the linearized, distributional forms 
of these identities in Fermi normal coordinates, we find A^^q' = nIj^ = A^]"^' = as well as 
doN^^^ = OoNI"^^ = and Mai = N^^^^d[oAf . As usual A^^^' and A^^^^ may be interpreted 
by their appearance in the body exterior metric and vector field. In the Einstein-Maxwell 
case discussed above one sees that the usual notions of mass M and charge q are related by 

use this to define q and M for a general theory with 
Maxwell gauge invariance, and rewrite the results covariantly to obtain 

= SvrM [ UauMx, z{T))dT (50) 

J-y 

^[A](l)a _ ^g^^ f (51) 

J 'y 

as well as 

Mu'WbUa = qu\2V[aA^^), (52) 

where q and M are constants. Thus for theories with the Maxwell gauge symmetry we have 
the usual point particle stress-energy and charge-current, along with the Lorentz force law. 
For theories without the gauge symmetry, however, the situation is more complicated. In 
this case we only have the single identity fH21) . and the Fermi coordinate calculation now 
gives 



n!^ = ^AtI^] 4°) (53) 
doNlfJ = -A^[^]^a[,4°) + ^A^^^doN^^^ (54) 
AiJa. = N^^^^d,,Af + IdoiNl^^A^^^) - ^^f ^^oA^^ ■ (55) 
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Since we no longer have NIq = Nl^ = Nl =0, the parameters q and M no longer suffice 
to characterize the body. In light of equation (l53l) . which requires A^^'^'j to point along Af\ 
it seems simplest to introduce a second charge q by A^'"^^, = lQTiqAf\ Then the results are 
rewritten covariantly as 

= 8^ / [MuaU, + AqP\^A,)Af) 5,{x, z{T))dT (56) 
E[A](i)a ^ ^ l^^^a ^ qpab^io)^ ^^^^^ z(r))dr (57) 

and 

(M - qiA'^A.y) u'^VaUb = 2{q - qA'u,)V ^,Aa]u'' 

+ P\ {2gA^V[,A,] + M^Ve (gA.M'^A^) - A.M'^Veg} (58) 
u'^WaM = -gA"M^(2V[aA5]) + wM^M^Vfeg, (59) 

where P'^^ = 6\ — ■u'^Mb projects orthogonally to u""^ and the superscript (0) on A^a'^ has 
been dropped in the last two equations for readability. Thus for theories without the gauge 
symmetry the usual point particle stress-energy and charge-current are not obtained (in 
that the distributional forms are not parallel to u"")^ the Lorentz force law is not the correct 
force law, and the particle is described by a time-dependent mass M as well and two charges 
q and g, neither of which has an evolution law (just as there was no law for the scalar 
charge in scalar-tensor theory). An example of a commonly studied theory to which these 
equations apply is Einstein-Proca theory. Of course, in most references Proca lagrangian is 
coupled to matter via an interaction term that by itself has the gauge symmetry, so that 
solutions to the Proca-matter system respect charge conservation (but charge conservation 
does not hold as an identity). If one restricts to such matter, presumably one would have 
g = and u"''Vaq = for bodies made of that matter, whence the usual particle equations of 
electromagnetism ( 150H52|) would be recovered. However, it is not clear that there is reason 
for matter to conserve charge in a vector theory beyond the fact that it is required in the 
most familiar vector theory. It should be emphasized that the Proca theory (and others 
without the gauge symmetry) admit far more general behavior in the motion of bodies than 
does ordinary electromagnetism. This general behavior would give the motion of any matter 
that did not conserve charge, as well as the motion of any non-matter objects (such as black 
holes or "geons" ) that might exist in the theory. 
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Note finally that the analysis of this section can be straightforwardly generalized to 
the case of multiple vector fields (or even multiple vector and scalar fields) in the manner 
discussed in section IIIII for scalar-tensor theory. In this case one simply obtains copies of 
equations fl5B]) and flFTI) for each field, and the right-hand-sides of flSSl) and fISU]) are copied 
for each field to form a sum. If any of the fields have Maxwell gauge-invariance, of course, 
the simpler terms from equations (15011521) may be used for that field. There may also be 
different gauge symmetries that provide different simplification. An important example is 
non-Abelian gauge theory, where the charges respect "gauge covariant" conservation. More 
precisely, if we label the set of vector fields by and their field equation operators by 
jT;ia ^ (^—g)^^/'^SS/SAl^ (using capital Latin indices for "gauge indices") then the gauge 
symmetry of the Lagrangian gives Va-E^" = Yli j f^"^^ E^°-A-^ as an identity, where f^"^^ 



are the structure constants as defined in 



13| (anti-symmetric on the first two indices). The 



usual Fermi normal coordinate calculation on this identity implies that the "hatted charge" 
vanishes for each body and gives an evolution law for the charges . The small body 
equations of motion are then 

Mu'VbUa = J2 /(2V[,AJ)m'' (60) 
u^VaM = J2 f^^q^AiAfu^'v!' (61) 

I,J,K 

«"V.g^ = 5^/"V^aV. (62) 

i,j 

The evolution is now fully deterministic on account of the extra symmetry of the Lagrangian. 
Note that if the structure constants are totally anti-symmetric, then M is constant and these 
reduce to "Wong's equations" [l^ . 



V. THE GENERAL CASE 



The procedure used for scalar and vector fields generalizes straightforwardly to higher- 
rank tensor fields. Here we motivate the scaling by analogy with the scalings used before. 
The scalings used for the metric, scalar, and vector fields were all such that the power 
of A cancelled powers of A resulting from the Jacobian of the coordinate transformation to 
scaled coordinates, allowing (|2]) and its scalar and vector analogs to hold. For a general-rank 
tensor field T'^i-"", , , I define the scaled version f^-'^^, , = x^'-'^T"^-''". , so that 

01... Dm' bi...bm bl...bm 
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the analog of (|2]) holds. Then the usual reasoning leads to the usual assumptions for each 
component of the tensor field, and one proceeds exactly as in the previous sections. I will 
summarize this procedure in the form of a proof of a theorem, below. It is straightforward 
to follow the steps to determine the force law for any theory of particular interest (although 
higher-rank fields are more seldom considered). 

Theorem. Let S be an action in four spacetime dimensions in the sense of appendix E of 
Wald J^y such that 1) the Lagrangian C depends differomorphsim-covariantly on the met- 
ric Qab and some set of tensor fields {ipi} (tensor indices suppressed), i.e., (f)*C{gab,'ipi) = 
C{(j)*gab,(f)*'ipi) is satisfied for diffeomorphisms 0; and 2) i?^^' = {—g)~'^^'^6S/6g'^^ and E^ = 
{—g)~^f'^5S/5il)i are (local) second-order differential operators on {gab^i^i}- Suppose there 
exists a one-parameter-family {gab{X),'ipi{X)} satisfying the analogs of the assumptions of 
section Wr^x Then, the worldline 7 (four-velocity u"" and four- acceleration a"") and a cer- 
tain function M defined on 7 satisfy equations of the form Ma'^ = and u°''VaM = F, 
where f^Ua = and both f and F are local tensor functions of u"" , a"", gab{S^)\-y, "0/(0)17, 



V(0/(O))|7, and certain tensor fields defined on 7. 

The proof is essentially to follow the steps of the previous sections. Since these steps 
are by now familiar, I will omit some details in the description here (allowing considerable 
savings on notation). Begin by varying S with respect to an infinitessimal diffeomorphism 
to derive the "Bianchi identity" . This results in an expression of the form 

V"^S = 5^ [ {E' @V^i),+ {VE^Q^j) J (63) 

where the notation {A © B)^ indicates a sum of terms, each of which consists of the tensor 
product of A and B contracted in some way to yield a dual vector. Note that the explicit 
form for an arbitrary-rank tensor field was worked out in 15|]. Now derive the effective 
point particle description. By assumption iii), the components of g^^J and are 0(l/r) (r 
near zero). By the assumption of second-order field equations and the analysis of appendix 
\K\ the components of eI^^^^^ and E^^^'^ as distributions are multiples of 6^{x^); take the 
coefficients to be Njf}{t) and N\t) (component indices suppressed), respectively. Now apply 



That is, adjoin the components of the ipi and ipi to those of the metric and rescaled metric (respectively) 
where they appear, and replace satisfaction of the Einstein equations with satisfaction of the field equations 
£;[9l = and = 0. 
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the linearized, distributional form of equation ( 163|) . By satisfaction of the background field 
equations (-E^fe^ = E^^^^ = 0), this identity takes the form 

V^Elfld) = ^ v4°))^ + (VE^^'^ V^f ) J , (64) 

^ u u 

where here and below V is the derivative operator associated with the background metric. 
We then have 

V^^i?] = E [(^' ® V^f 1,)^ + (VAT^ I,) J . (65) 

^ u u 

Adopting Fermi normal coordinates, one has 

5(3) (f ) [-doNoo + aiNio - -F] + (f ) [iV,o - 6^.] = (66) 
5(3)(f) l-doNo, + a^N,, + a,Noo - H^] + d,6^''\x) [N,, - /C,,] = 0, (67) 

where J-", Qi, Hi, and ICij are local functions of A^^, ipf^lf, dipf^\^, ai. Naming A^oo = it 
now follows from fl66|) and fl67|) that 

doM = aiQi + (68) 
Moj = 9o^j - ajJCij + T^i. (69) 

These are the small-body equations of motion for the theory (expressed in Fermi normal 
coordinates). It is straightforward to determine the functions J-", Qi, Hi, and ICij by direct 
calculation for the theory in question (or even the general case; however, the expression 
is not simple). It is also straightforward to relate the parameters M and A^^ to the body 
exterior configuration {(^^p* , ^1°''} for each particular theory (as done in previous sections), 
enabling their calculation for any particular body. The covariant translation of equations 
fl68|) and fl69|) proves the theorem. 

The theorem establishes that a simplified description of motion (in the form of a second- 



order equation for 7)[^ is obtained via buffer- zone dynamics in a very large class of theories. 
It is interesting to speculate on the extent to which the class of theories could be enlarged. 
The requirement of second-order field equations seems easiest to relax, since higher-order 
equations would lead only to higher-order delta functions (i.e., derivatives of delta functions) 
appearing in the effective stress-energy, whence our calculations could proceed straightfor- 
wardly. However, it is far from obvious that our assumptions on one-parameter-families 



In some theories the field equations may enforce M = in which case the equation for 7 may be lower 
than second-order or even trivial. 
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would remain reasonable in the context of higher-order theories, whose solutions may have 
very different propertieso The requirement of a diffeomorphism-covariant Lagrangian ap- 
pears difficult if not impossible to relax, since the Bianchi identity plays an essential role 
in determining the motion. At least in our approach, the diffeomorphism-symmetrry and 
Lagrangian formulation are key to obtaining a description of motion from the buffer-zone 
field equations alone. 

VI. SUMMARY 

I have treated the motion of small bodies in classical field theory via the approach of 
J]. The search for one-parameter-families of solutions representing the exterior field of a 
shrinking body led precisely to the physical assumption of a "buffer zone" — a region far 
enough from the body that its field can be approximated in a multipole series, but close 
enough to the body that the field of the external universe can be approximated in an ordinary 
Taylor series. No assumptions about the body interior are made. In the case of second-order 
metric-based theories following from a diffeomorphism-covariant Lagrangian, I derived the 
force law for scalar and vector fields, and showed that the method works in the general-rank 
case. This provides a rigorous derivation of the small-body force law in many classical field 
theories commonly considered, and shows that field dynamics outside a body determines its 
motion in a very general class of theories. 

ACKNOWLEDGEMENTS 

I wish to thank Wayne Hu, Fabian Schmidt, Michael Seifert and especially Bob Wald for 
very helpful conversations. This research was supported in part by NSF grants PHY04-56619 
and PHY08-54807 to the University of Chicago. 

Appendix A: Delta-function Calculation 

Consider a linear, second-order differential operator L that takes tensors of rank 

n into tensors of rank m. The adjoint is the linear map from tensors of rank 



^•^ Note, however, that many higher-order theories (such as those whose Lagrangian is a function of the Ricci 
scalar) admit second-order formulations, so that the current analysis applies. 
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m to tensors of rank n defined by 16 1 

for arbitrary ip and 0. The vector s" is a multilinear function of 0«i-any^^bi-bn^ ^bi...6„y^0ai...a„ 
and (j)°-i---°-"^bi...bn _ j^^^ promote L to an operator on distributions in the standard way. 
That is, we define a distributional operator L by 

{U)[f] ^ I d'xV^{L^fr-''-<jya,,„a^ (A2) 

for smooth test tensors /^i---^'" of compact support. Now suppose that the metric and 0-field 
have the following expressions in coordinates 

g^.u = Vf^u + 0(r) (A3) 
r 

and further that ^'^i---^'" solves the equation at r > 0, i.e., that we have 

(^L(f)fi-^m. _ 0, r > 0. (A5) 

Here r, 6*, are related to in the usual way. The distribution [LcpY^'"'^"' may be computed 

by 

{Urn = hm / d'xV^{L^fr-^"<j>,,...,^ (A6) 



e-S>0 



hm / d'x^ + VpS"} (A7) 



^ ^ I — I > -i I " "I I .f 1 1 I — I — » / 

e>0 



= lim / sin 9d9d(t)dt UpS^, (A8) 

where Up is the unit normal to the unit two-sphere in Euclidean three-space, and is 
evaluated on and /. The first line follows from (lA2p and the smoothess of f^^^■■■^^", The 
second line follows from ( lAll) . The third line follows from ( ]A5I) and integration by parts, 
where the boundary term vanishes by the compact support of J/^i - a*". The volume element 
on the surface r = e has been replaced with the Minkowski volume element on account of 
the limit e — )■ and the metric form flA3l) . (Note that if different coordinates were chosen 
such that the metric is not Minkowski at r = 0, the volume element would take a different 
form, but it is easy to check that the analysis would still hold.) Now, equation flA4p together 
with the properties of imply that s'' may be written 

s' = \D'^''^-''-{t,eA)U...uJt,r,e,<j>) + 0{l/r) (A9) 



26 



for some Df^^-'^"'. Then we have 

(L0)[/] = hm / sin 9d9d(f>dtnpDP''' -""^f^, ^Jr=e (AlO) 
dtU..,Jr=Q I sinOded^ ripDP'''-''-^ (All) 



= j dtN^^-^-{t)U,....^{t, x' = 0) (A12) 

where the second step follows from the smoothness of f^-'^rn^ a^^^j ^j^g i^^^ g^gp simply defines 
the result of the angular integral to be N^'^-^"'{t). Since the test function is evaluated at 
X* = 0, this expression shows that distribution {Lcj))"^-'''^ is proportional to the spatial delta 
function 5*^^)(a;); i.e., we have the desired result 

(^LfPyi-'^m ^ Ar'^i-'^-(i)(5{3)(f). (A13) 

Note that an explicit formula for N'^^-'^"'{t) in terms of C^^...ju„(i, ^, 0) can be determined 
for any particular differential operator L by following the steps of this computation explic- 
itly. This formula will in general involve angular integrals of C/i^...^^ and its first angular 
derivatives. 



Appendix B: Scale Invariance 

In this appendix I provide a definition of scale invariance and compare with other defini- 
tions commonly given. Tensor indices are suppressed throughout this section. It is conve- 
nient to restrict to theories that follow from a Lagrangian, 

S(g,^l;j)^ j d'xC{gM, (Bl) 

although the metric need not be dynamical and there may exist other "background struc- 
ture" . (That is, the Lagrangian is not required to depend diffeomorphism-covariantly on g 
and The theory is scale invariant if the action scales homogeneously under a scaling of 
the metric and other fields; more precisely, if exist numbers {-P/} and n such that 

S{\^g,\'''i;:)^\^S{gM, (B2) 

for numbers A. This property implies that the equations of motion for the rescaled fields are 
identical to the equations of motion for the original fields, so that there is no preferred scale 
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for lengths (as measured by g) or field values (0) in the theory. The choice of for the 
scaling of g is conventional, and gives A the interpretation of a length (since the line element 
scales as A^). Note that our definition of ^ = X^^g in the body of the paper is consistent 
with this convention of "(7 — X^g" . 

An alternative definition of scale-invariance is often given in the context of theories spec- 
ified by partial differential equations in coordinates (without the introduction of a metric). 
In this case one rescales the coordinates Xx^ and asks if a rescaling of the fields 

can restore the original partial differential equations. Such theories can often be rewritten 
diffeomorphism-covariantly via the introduction of a fiat metric gab- Changing the coordi- 
nates of the original partial differential equation then corresponds to applying the diffeomor- 
phism a associated with coordinate rescaling to all fields except the metric. That is, to check 
if a diffeomorphism-covariant equation E[g,il)i\ = is scale-invariant according to the non- 
tensorial definition, one asks if a^E[g,%lJi\ = X'^E[g, X^'a^ipi]. Applying cr^^ to both sides 
and using the diffeomorphism-covariance of E, this becomes E[g,ilJi] = X"'E[a~\g, X^'ipj]. 
But since g is fiat, cr~\g = X^g, and this reduces to the definition involving rescaling the 
metric (at the level of the field equations). This is why one says that rescaling the metric is 
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a curved-spacetime generalization of rescaling the coordinates 

In the context of quantum field theory a second alternative definition of classical scale- 
invariance is often given by requiring that the action be left invariant under the rescalings, 
rather than just scale homogeneously 



18| . This corresponds to our definition with the 
additional demand that n = 0. Since the overall scaling of the action (i.e., the value of n) 
has no effect on the classical equations of motion, this notion in fact removes as "classically 
scale-invariant" many classical theories with no preferred scale. For example, vacuum general 
relativity is scale-invariant in our sense, having no preferred length scale; however, vacuum 
general relativity is not classically scale invariant according to the definition used in quantum 
field theory. Another example is massless 0" (n 7^ 2) in flat spacetime (scale-invariant in 
our sense); only 0^ is classically scale- invariant in the sense used in quantum field theory. 

The difference between the definitions can be further elucidated with reference to the 
well-known fact that a scale-invariant action contains only coupling constants that are di- 
mensionless in particle physics units (c=l, ^=1, Gt^I). Adopting the viewpoint that 
masses and lengths are fundamentally different (but that time intervals and lengths are not), 
it is convenient to work in special relativity units (c = l,h 7^ 1,G 7^ 1), where the statement 
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of dimensionlessness in particle physics units becomes the property of having equal mass 
and length dimension, so that in particular one cannot construct a length using h (which 
has dimensions of mass times length). For example, the constant A of A(f)^ has dimensions 
of mass times length so that no length can be constructed with h. On the other hand, the 
constant G of general relativity has dimensions of mass over length, so that h can be used to 
construct a length (called the Planck length) . While our notion of classical scale-invariance 
implies that no lengths can be constructed from the coupling constants alone, the quantum 
field theory notion of classical scale-invariance places the further restriction that no lengths 
can be constructed even when one is allowed to use h in addition. 
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